National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 


Best 

Available 

Copy 


UN  CLASSIFIED 

SECURITY  CL  ASSI  FIC  ATIOH  OF  This  PACE  (irh*n  Dele  Entered) 

REPORT  DOCUMENTATION  PAGE 


)//1C'Q2-<2£:-/L 

READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM  . 


I REPORT  NUMBER 


'2.  GOVT  ACCESSION  NO.i  >•  RECIPIENT'S  CATALOG  NUMBER 


« TITLE  (end  Submit) 

UNSTEADY  HYDRODYNAMICS  OF  A 
BODY  OF  REVOLUTION  WITH 
FAIRWATER  AND  RUDDER 

4 AUTHOHflJ 


Douglas  P,  G las  son 

’9  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Massachusetts  Institute  of  Technology 
Department  of  Ocean  Engineering 
Cambridge,  Mass.  02139 

II.  CONTROLLING  OFFICE  NAME  ANO  AODRESS 


5.  TYPE  OF  REPORT  & PERIOD  COVERED 

Final  Report 

Oct.  71-March  74 

6,  PERFORMING  ORG.  REPORT  NUMBER 


a CONTRACT  OR  grant  nuubero; 

N00014-67-A-02  0 4-0058 


10.  PROGRAM  ELEMEN  *.  PROJECT,  TASK 
AREA  ft  WO'*’*  UNIT  NUMBERS 


12.  REPORT  DATE 

March  74 

13.  NUMBER  OF  PAGE^ 


14  MONITORING  AGENCY  NAME  A ADDRESS*”//  d/l/arent  Irvn  Controlling  Ofllce)  15.  SECURITY  CLASS,  (ol  thla  report) 

UNCLASSIFIED 


15a.  DECLASSIFICATION/DOWNGRADING 

schedule 


i 16.  DISTRIBUTION  STATEMENT  (of  thla  Report) 


This  document  has  been  approved  for  public  release  and  sale;  its 
distribution  is  unlimited 


17.  DISTRIBUTION  STATEMENT  (ol  tha  abstract  mntarad  In  Block  20,  It  dl  frrent  trom  Raport) 


18  SUPPLEMENT  ARY  NOTES 


NATIONAL  TECHNICAL 

information  service 

US  *DA'*Ti#.r.l  rf  f 
'-I-F  • 3 *...!<!  VA  ? ? 1 b » 

19  KEY  WORDS  (Continue  on  rmvaraa  aide  It  neceaeatry  and  Identify  by  block  number) 

UNSTEADY  LI  FTING  S UR  FA  CE 
SUBMARINE  CONTROL 
SLENDER  BODY 

20  ABSTRACT  (Continue  on  revetae  aide  If  nateaaary  and  Identify  by  block  mmbor) 

Potential  flow  models  arc  developed  for  a submerged  body  of 
revolution  with  fin  and  rudder  appendages. 

Forces  and  moments  on  the  lifting  surfaces  and  hull  have  been 
predicted  at  a steady  angle  of  attack.  The  procedure  is  extended  to 
the  time  dependent  angle  of  attack  case.  Experimental,  analytical  and 
numerical  approaches  are  described. 


DD  i j an 1473  EDITION  OF  I NOV  65  IS  OBSOLETE  * 

I , 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fHTion  Del • Entered) 


-JrV707'77^^Tf-T?WV»<l 


MASSACHUSFTTS  INSTITUTE  OF  TFCHNO LOG Y 
Department  of  Ocean  engineering 


Report  74-7 

UNSTFADY  HYDRODYNAMICS  OF  A 
BODY  OF  RF VOLUTION  WITH 
FAIR  WATFR  AND  RUDDFR 


Douglas  P.  Glasson 
March  1974 


Approved  for  public  release; 
distribution  unlimited 


This  research  was  carried  out  under  the  N ,/al  Ship  Systems 
Command  General  Hydromechanics  Research  Program  Sub- 
project SR  009  01  01,  administered  by  the  Naval  Ship  Research 
and  Development  Center.  Prepared  under  O.  N.R.  contract 
N000 14-67 -A -0204-0058.  M.I.T.  DSR  73836 


II 


Abstract 


Potential  flow  models  are  developed  for  a submerged  body  of 
revolution  with  fin  and  rudder  appendages. 

Forces  and  moments  on  the  lifting  surfaces  and  hull  have  been 
predicted  at  a steady  angle  of  attack.  The  procedure  is  extended  to 
the  time  dependent  angle  of  attack  case.  Experimental,  analytical  and 
numerical  approaches  are  described. 
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I nt  roduction 


The  primary  aim  of  the  work  presented  in  this  report  is  to  add 
insight  i n t <j  dynamic  interaction  effects  experienced  by  high  speed 
submarines.  The  dramatic  advances  in  submarine  design  since  the 
second  World  War  have  been  in  the  direction  of  higher  speed,  larger 
size  and  endurance,  but  depth  capability,  expressed  in  boat  lengths 
or  seconds  has  not  increased  proportionally.  As  a result,  submarine  control 
design,  particularly  for  vertical  plane  motions  and  c oupling  between 
maneuvering  control  in  the  horizontal  plane  and  response  in  the  vertical 
plane  has  become  of  great  importance.  Its  importance  will  urdoubtably 
increase  as  submarine  technology  develops. 

As  a result  of  the  demands  on  submarine  control  capability, 
unc1'*  rstanding  of  the  hydrodynamic  forces  acting  on  the  boat  has  assumed 
greater  importance  and  effects  which  were  considered  negligible  or  merely 
aggravating  at  low  speed  may  be  crucial  at  high  speed. 

The  usual  approach  to  solving  problems  in  marine  vehicle  dynamics  is  to 
express  the  forces  acting  ori  the  vehicle  due  to  its  motion  as  functions  of 
its  velocity  and  time  derivatives  of  velocity.  Such  an  approach  is  extremely 
useful  since  it  lends  itself  conveniently  to  linearization  about  discrete 
operating  points  of  interest  and  usually  higher  order  derivative  terms  are 
of  diminishing  importance. 

There  are,  however,  disadvantages  to  this  general  approach. 

When  the  vehicle  alters  its  pnviro  nment  in  such  a manner  that  the  forces 
and  moments  acting  upon  it  are  dependent  on  its  motion  history  over  a 
significant  time  span,  higher  time  derivatives  become  more  important  and 
the  coefficients  of  these  time  derivatives  become  difficult  to  predict  by 
theoretical  or  experimental  techniques. 

An  example  of  such  an  effect  on  the  vehicle  environment  is  the 
trailing  vortex  system  of  a submarine  sail.  The  relative  location  and  strength 
of  this  vortex  system  depends  on  the  circulation  history  of  th  sail,  the 


motion  of  the  vortex  sheet,  anti  the  trajectory  and  orientation  of 

the  vehicle.  The  vortex  sheet  has  an  effect  on  the  hull  and  rudder  forces 

which  is  strongly  dependent  on  the  relative  location  of  the  trailing  sheet. 

Theiefore  the  forces  on  the  vehicle  are  dependent  not  only  on  its  present 
velocity  anti  acceleration,  but  on  its  motion  history.  To  express  this 
effect  within  the  traditional  framework  demands  enough  higher  time 
derivatives  of  the  motion  to  recreate  the  vehicle  history  between  the  time 
the  circulation  formed  on  the  sail  to  when  it  ceases  to  effect  the  rudder. 

There  is  of  course  no  guarantee  that  higher  derivatives  will  be  of  monotonica lly 
diminishing  importance. 

The  approach  taken  in  this  work  is  to  consider  the  hydrodynamic 
forces  acting  on  the  vehicle  directly  as  a function  of  its  motion  history.  The 
interactions  between  vortex  sheets,  hull,  and  rudder  are  derived  directly 
from  the  geometry.  This  implies  that  for  simulation  work  the  motion  history 
of  vehicle  and  vortex  sheets  must  be  stored  as  the  alternative  to  employing 
a large  number  of  time  derivatives. 

The  thrust  of  this  effort  is  to  develop  techniques  for  evaluating  hydrodynamic 
restoring  forces  in  a history  dependent  situation.  Various  techniques  have 
been  employed,  but  the  most  promising  for  further  pursuit  seems  to  oe 
a slender  body  approximation  for  the  hull  with  singularity  distribution 
models  for  the  lifting  surfaces.  In  both  the  slender  body  and  lifting  surface 
calculations,  rough  approximations  and  widely  spaced  grids  are  used. 

The  purpose  of  the  authors  was  to  develop  techniques  and  evaluate 
the  importance  of  effects  not  normally  taken  into  account.  Therefore  little 
effort  has  been  invested  in  refining  numerical  results.  The  difference  between 
results  including  and  not  including  interactions  is  the  point  of  interest.  If 
the  techniques  investigated  are  to  be  used  for  simulation  or  prediction  of 
submarine  motions,  the  numerical  procedures  described  require  considerable 
refinement . 


J 


The  total  effort  divides  into  several  logical  parts.  First  an 
experimental  study  was  conducted  in  the  MIT  variable  pressure  water 
tunnel  of  a sma'l  submarine  model  with  various  fin  configurations.  Forces 
and  moments  were  measured  at  various  angles  of  attack.  Of  greatest 
importance  in  the  experimental  study,  liowever,  werethe  visual  results. 

The  facility  used  has  the  great  advantage  of  plexiglass  walls  and  control 
of  ambient  pressure.  By  lowering  the  tunnel  pressure,  cavitation  bubbles 
can  be  induced  in  the  trailing  vortex  system  from  the  sail  so  that  its 
trajectory  may  be  observed.  In  addition  tufts  may  be  attached  to  the 
hull  and  fins  to  observe  flow  patterns  around  the  model. 

These  studies  are  reported  by  Luckard  (1)  under  separate  cover. 

An  analytical  approach  to  the  interaction  problem  was  pursued  by 
Newman  and  Rodriugez  making  slender  body  and  low-aspect  ratio  assumptions 
and  linearizing  the  problem.  This  approach  is  described  by  Newman  and 
Wu  i 2). 

The  bulk  of  the  work  was  devoted  to  numerical  approaches  to 
representation  of  the  hydrodynamics  problems.  The  approach  assumes  that 
the  hull  is  slender  and  a body  of  revolution  but  includes  finite  aspect  ratio 
lifting  surfaces.  Thr  boundary  value  problems  are  approximately  solved  with 
d;screte  singularity  distribution  representations  of  hull,  lifting  surfaces  and 
rotational  wake. 

Experimental  Approach 

The  single  greatest  question 
where  the  trailing  vortex  wake  from 
approximations  to  the  hull  forces  as 
necessary  to  assume  that  the  vortex 
angle  of  attack.  For  small  sideslip 
the  sail  is  assumed  to  remain  in  the 


at  the  initiation  of  this  work  was 
the  sail  went.  To  obtain  linear 
a function  of  sideslip  angle  it  is 
sheet  trajectory  is  not  dependent  on 
angles  the  trailing  vortex  sheet  from 
plane  of  the  rudder  axis  and  longitudinal 


cent  e rl  m.- . In  this  Inration  the-  sa  1 wake  has  maximum  effect  on  rurlder 
lift.  If  in  practice  the  sail  wake  does  not  pass  close  to  the  rudder,  this 
effect  will  be  appreciably  reduced.  The  path  of  this  trailing  vortex  sheet 
is  infb  enced  by  its  own  induce*.,  velocities  as  well  as  the  free  stream 
and  the  vehicle  boundaries.  Approximati ons  may  be  made  on  the  basis  of 
lifting  line  theory  and  singularity  distributions  fer  the1  trajectory,  but  the 
tendency  of  the  sheet  to  roll  up  makes  such  computation  questionable. 

For  this  reason  photography  of  actual  vortex  sheet  paths  at  steady 
angles  of  attack  was  valuable  input  to  the  modeling. 

The  model  used  for  the  experimental  and  numerical  work  is  two 
foot  long  submarine  like  ' body  of  revolution  equipped  with  removable 
fairwater  and  upper  and  lower  fins  for  testing  in  the  MIT  variable  pressure 
water  tunnel.  Perhaps  the  most  important  result  of  this  experimental 
and  analytical  effort  is  that,  at  small  angles  of  attack,  the  lift  on  the  rudder 
behind  the  sail  is  -oversed  from  that  predicted  if  the  wake  of  the  fairwater 
is  ignored.  This  is  of  importance  primarily  to  directional  stability  pre- 
dictions of  submarine  type  hulls.  The  effect  of  the  rudder  is  destabilizing 
rather  than  stabilizing  due  to  its  interaction  with  the  fairwater  wake.  The 
most  recent  work  on  the  unsteady  hydrodynamics  problem  comprises  the 
technical  content  of  this  report  since  the  experimental  and  steady  numerical 
approaches  are  included  in  the  report  by  Luckard  (1).  Comparisons  of 
results  for  the  steady  case  with  the  Newman  and  Wu  (2)  analytical 
results  are  also  presented  by  Luckard  (1). 

The  part  of  this  project  that  is  probably  of  most  general  interest  is 
the  approach  taken  to  numerical  solutions  of  the  forces  on  fairwater,  hull, 
and  rudder  due  to  a sudden  change  in  angle  of  attack. 
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I!.  The  Response  of  a Submarine  Sail  to  a Sudden  Change  in 

Slides  lip  \ ngl  e 


A.  Unsteady  Lifting  Surface  Theory 


The  transient  buildup  of  lift  and  moment  on  an  initially  unloaded 


'ifting  surface  that  has  undergone  a sudden  change  in  angle  of  attack 

is  a direct  result  of  the  physical  fact  that  a finite  time  is  required  for  the 


trailing  vortex  sheet  of  the  lifting  surface  to  attain  ii  s steady  state 


configuration.  At  the  instant  that  the  surface's  orientation  to  the  flow  is 


changed,  all  vorticity  is  confined  to  the  surface  itself.  As  the  surface 


moves  forward  in  its  new  orientation,  vorticity  is  shed  from  the  surface 
into  the  external  flow.  Eventually,  the  lifting  surface  reaches  a steady 


state  condition  in  which  the  wake  has  attained  its  familiar  "trailing  vortex 


sheet  "'  st  rucui re , except  for  regions  of  the  wake  that  can  be  considered  to  be 


an  infinite  distance  from  the  lifting  surface. 


Using  a distribution  of  dipoles  to  represent  the  lifting  surface 
and  its  wake  at  a given  instant  of  time,  the  component  of  velocity  is  given  by 
the  following  expression:  (Derivation  is  in  Appendix  A) 
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where  X is  the  local  dipole  sheet  strength,  its  second  mixed  derivative 


expressi  lg  the  local  strength  of  the  rel  ated  point-horseshoe  vortex  sheet. 


Applying  the  boundary  condition  that  there  is  no  flow  normal  to  the  lifting 


se  face  boundary  at  E»Oone  obtains: 
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The  integral  over  the  wake  has  been  written  in  this  manner  to 
allow  the  option  of  placing  the  wake  in  a plane  other  than  /.  0.  Equation 

(II.  2)  is  the  mathematical  statement  of  the  relationship  between  the  known 
boundary  conditions  for  the  lifting  surface  and  the  unknown  singularity 
distribution  representing  the  lifting  surface  and  its  wake. 

The  irrotationality  condition  on  the  flow  external  to  the  lifting  surface 
and  its  wake  generates  an  important  relationship  between  the  distribution  of 
vorticity  in  the  wake  and  the  history  of  the  total  vorticity  bound  to  the  lifting 
surface.  Consider  a two  dimensional  section  of  the  wing-wake  flow.  In  order 
for  the  value  of  £■$  (the  local  dipole  sheet  strength)  to  be  single  valued  at  the 
trailing  edge  of  the  lifting  surface: 


Since 


is  the  total  vorticity  strength  enclosed  by  the 


Ca 

contour  and,  hence,  the  total  bound  vorticity  on  the  wing  at  the  moment 
of  interest,  this  relationship  states  that  the  total  amount  of  vorticity  shed 
into  the  wake  equals  the  negative  of  the  present  amount  of  vorticity  bound 
to  the  lifting  surface..  Consider  this  process  of  bound  vorticity  increase  and 
shedding  of  vorticity  into  the  wake: 
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If,  ai  any  section,  the  total  amount  of  bound  vorticity  increases 
by  ar„  in  dt  amount  of  time,  then  an  amount  of  circulation  equal  to 
must  have  been  shed  into  the  wake  such  that: 
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Substituting  into  (II. 2),  the  boundary  value  integral  equation  becomes: 
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Hence,  the  integral  over  the  wake  has  been  reduced  to  a convolution 
integral  involving  (known)  previous  values  of  the  total  bound  circulation 
at  each  section  of  the  lifting  surface. 

In  order  to  calculate  the  forces  and  moments  on  the  lifting  surface 
as  a function  of  distance  traveled,  it  :s  necessary  to  consider  the  time- 
dependent  form  of  the  Bernoulli  equation. 
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whe  re  CD  and  C£>  denote  points  on  a streamline 

Applying  this  equation  to  points  on  the  upper  and  lower  surfaces  of  a 
lifting  surface  section  with  a steady,  uniform  free  stream  one  obtains: 
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This  expression  can  also  be  written: 
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is  then  integrated  in  the  usual  manner  to  yield  the  total  force 
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and  moment  on  the  lifting  surface. 


Although  a complete  solution  of  the-  problem  of  the  transient  response 
of  a lifting  surface  requires  a solution  of  equation  (II.  6),  some  important 
trends  ca.i  oe  deduced  by  considering  the  startin'’  instant  and  steady 
state  solutions  of  equation  (II.  fl.  For  a flat,  rectangular  wing  of  unit 
half  chord  with  a planar  wake,  equation  (II  , ( I reduces  to- 
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It  has  been  shown  by  Wagner  ( 1 that  in  the  two  dimensional  {&■*  •*) 

limits  of  the  above  equations  the  starting  problem  is  satisfied  by  a 
vorticitv  mode  of  tr  functional  form 
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I’he  rosuF  for  the'  steady  two-dirr  nsional  case  from  classical 


thin  airfoil  theory  s- 
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The  sectional  1 i f t for  the  starting  problem  is,  Ref.  (1): 


CL  - Tio/ 

h\r«k 


ft.  IT) 


and  f cr  the  steady  problem: 
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For  finite  aspect  ratio  wings  the  presence  of  the  trai’mg  vortex 


wake  and  its  associated  downwash  causes  a decrease  in  the  ratio  of  the 


final  and  initial  lifts.  As  the  aspect  ratio  of  a finite  wing  is  decreased, 
the  waki  makes  an  increasingly  dominant  contribution  to  the  boundary 
condition  equation  (II.  6),  until  in  the  limit  of  zero  aspect  ratio,  the  boundary 
condition  is  assumed  to  be  satisfied  exclusively  by  the  trailing  vortex  sheet. 

It  seems  reasonable  that  for  sufficiently  low  aspect  ratio  the  starting 
lift  may  exceed  the  steady  state  lift  due  to  the  predominance  of  the  trailing 
vortex  sheet  downwash  in  the  steady  case. 

A numerical  study  was  carried  out  to  compare  the  starting  and 
steady  lifts  of  isolated  rectangular  wings  of  various  aspect  ratioes. 

The  .ntegral  equation  (II.  9)  was  solved  for  the  starting  problem  assuming  a loading 
expansion  of  the  form: 
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1 e.  the  chordwise  mode  shape  for  the  two-dimensional  case  multiplied 
by  some  spanwise  distribution  function'.  In  the  steady  case,  the  loading 
t xpansion  assumed  was  of  the  form: 


ILe  t,  Z S;l*j)  (J3.t)v 
Bx  r-i-i  - .Jj,  ' 

C*lx)«  "1l+v  > u.«*5 
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and  the  Si' are  given  by: 


© 3 SlWJ 


In  all  cases,  the  unknown  coefficients  of  the  assumed  loading  modes 
were  obtained  by  applying  equation  (II.  f)  at  16  evenly  distributed  control 
points  and  finding  a best- lea st- squa re s fit  of  the  J i to  the  (over  determined! 
set  of  linear  equations  that  resulted.  The  starting  and  steady  forces  were 
then  calculated  by  integrating  the  sectional  forces  given  by  equations  (11.13) 
and  (11.14)  across  the  span. 

The  resulting  values  of  Cj  steady/Cj  start  are  plotted  below'.  The 
results  show  that  the  starting  and  steady  lifts  are  approximately  equal  at  an 
aspect  ratio  of  2.  At  an  aspect  ratio  of  1.  (a  typical  value  for  the  aspect 
ratio  of  a modern  submarine  sail)  the  starting  lift  exceeds  Ihe  final  lift  by 
roughly  2 3%.  Hence,  a significant  quantitative  difference  between  the  forces 
calculated  in  the  transient  (step  response)  of  a low  (“Si)  aspect  ratio  wing 
and  what  would  be  predicted  in  a psuedo- steady  analysis  is  already  apparent. 


R.  Indicia)  Response  of  a Submarine  Sail: 


1)  The  Image  System 

Since  the  sail  is  attached  to  what  often  can  be  considered  a locally 
i yiindrical  surface  (the  submarine  hull),  an  image  system  is  required  to 
approximately  satisfy  the  boundary  condition  of  flow  tangency  to  the  hull. 
The  theoretical  basis  for  the  image  system  used  is  given  in  Milne  - 
Thompson  ( 5 ) and  is  i eviewed  by  Luckard  ( 1 ) for  the  particular  case  of 
a hull-sail  combination.  Using  Luckard' s results  for  the  calculation  of 
the  equivalent  span  of  the  sail,  one  obtaii  s: 


where  RSI.  and  RMS  are  given  in  terms  of  half  chords  of  the  sail.  In 
addition  to  determining  the  equivalent  span,  the  presence  of  the  hull  causes 
an  increase  in  the  in-flow  velocity  to  the  sail  given  by  ( 1 ) ; 

A?  ~ \JooB 


Hence,  the  local  angle  of  attack  seen  at  a control  point  on  the  sail 

i s gi ven  by : 
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) Th<*  Spanwiso  Loading  Modes: 


Fhe  tangency  condition  (on  the  hull  surface)  also  requires  that 


the  sheet  of  trailing  vorticity  have  zero  strength  at  the  hull  surface, 


i.e.  the  discontinuity  in  the  y perturbation  velocity  across  the  trailing 


vortex  sheet  must  vanish  at  the  hull-sail  junction.  This  reauires  that 


the  slopes  of  the  spa.iv.ise  loading  functions  assumed  to  solve  equation  (II.  6) 
be  zero  at  the  hull-sail  junction.  To  satisfy  this  requirement,  spanwise 
loading  modes  of  the  forms  given  below  were  chosen: 


. vw  lx-1  • ; vh‘-V-1 

(II. M) 

■ Mi'i]  ' 

uhc  re:  . r £s,V  J OZ  ^ * f ^ [f) 


where: 


3)  The  Chordwise  Loading  Modes 


In  addition  to  the  chordwise  loading  modes  given  by  (11.17)  a 
mode  is  requ.red  that  can  provide  for  a non-z.cro  value  of  vorticity  at 
the  t railing  edg<  of  the  sail . Since  the  starting  mode  given  by  (II . 11)  yields 


a nearl  constant  chordwise  downwash  for  the  aspect  ratio  range  of 


i it o rp st , a special  mode  was  constructed  having  the  following 


functional  form: 


{1?*)  i+%  i i"  iur 


Ox) 


131.20) 


*-%\i 

i*-*vd 


H B . ofa^5 
2. 


The  first  term  represents  the  chordwise  vorticity  distribution  given  by 


(II.  11)  where  the  singular  trailing  edge  portion  of  the  distribution  has  been 


moved  past  the  trailing  edge  of  the  chord  by  replacing  X with 


( K-% 

\~OT. 


The  second  term  is  subtracted  in  order  to  make  the  net  circulation  of 


the  special  mode  equal  to  zero. 

Since  equation  (II. 4)  relates  the  trailing  edge  value  of  vorticity 


to  the  total  net  vorticity  bound  to  the  lifting  surface,  it  is  convenient  to 


normalize  the  trailing  edge  value  of  vorticity  for  this  chordwise  mode  by 


dividing  through  by  C^fl] 


I - (—  T 

r.x  (li  \a  —1™  _ j* 

Clx)=  \i-%)]  1-  irU-Hj 


(n.zj) 


T — — • 


4)  The  V ake 

The  distribution  of  vorticity  in  the  wake  at  any  point  in  time  was 
assumed  to  have  a y variation  that  eould  be  described  by  the  spanwise 
modes  (11.191  and  a piecewise  linear  x variation.  The  piecewise  linear 
variation  in  x was  chosen  since  this  reduces  the  convolution  integral 
representing  the  downwash  of  the  wake  in  equation  (II.  6)  to  a simple 
matrix  multiplication  of  a dovnwash  influence  matrix  by  the  previous 
values  of  the  total  bound  vorticity  on  the  lifting  surface.  A piece  wise 
linear  distribution  also  has  the  advantage  that  the  x integration  of  the 
second  term  of  equation  (II.  6)  can  be  carried  out  analytically. 

The  trailing  edge  of  the  wake  for  the  indicial  lift  case  has  a 
square-root  singularity  of  vorticity  strength  in  x.  Since  this  distribution 
can  be  comparatively  troublesome  to  integrate  numerically  or  otherwise, 
this  portion  of  the  wake  was  modeled  by  a piecew-ise  linear  portion  plus  an 
impulse  of  vorticity  at  the  trailing  edge.  The  relative  values  of  the 
piecewise  linear  portion  and  the  impulse  were  adjusted  to  make  the  areas 
centroids  of  thetwo  distributions  equal. 


^rp*rr--yj-’~}- ‘^TVT 7-J^>Tf£7^*^M ./ I l**u  iv«y  r TT*.  >y *5,',P”  M',*’i1^ r~- WT ’ .'V  TV7*  rr'.^-tj’rflj 

..  ..  ■«-3 
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5)  Step  by  Step  Solution: 

At  this  point,  equation  (II.  6)  can  be  written  as: 


^lit)  ls-0 


^(Xojtjojb)  “ C>w*(i  ) ^WS^lXo^o)  f D^JlK.iXojIja^T'  ^ 

a 

+ ^ -31^(0 J (u.a) 


where  fi  is  the  local  sideslip  angle;  bbtlS^,*,  (,  Xo  jl^o  ^ is  the  side  wash 
at  control  point  ( ( \^„  ) on  the  sail  due  to  a vorticity  cross  mode  m 

chordwise  - n spanwise.  is  the  coefficient  of  crossmode  mn  . 

r ) is  the  sidewash  at  control  point  Ix.j  on  the  sail  due 

to  a spanwise  strip  of  the  wake  vorticity  distribution  located  a distance  X 
from  the  trailing  edge  of  the  sail.  DUW.kOc*,^)  is  the  sidewash  due  to 
the  special  chordwise  mode  given  by  (11.21).  This  equation  is  to  be  solved  for 
a discrete  set  of  distances  traveled  s. 

Since  the  distribution  of  vorticity  in  the  wake  of  the  sail  was  assumed  to  be 
piecewise  linear,  relationship  (11.4)  must  be  finitir-.cd . Consider  the  balance 


of  wake  and  bound  vorticity  at  a section  of  the  sail  during  a finite  interval 


H 


of  1 rme  — with  the  assumption  that  the  vorticity  of  the  wake  is  linear 

•i  H.  Relationship  (II.  3)  requires  that: 


Uo 

between  x 1 and  x 


I if  i -r.o  3.t= 


Ac  = - ^ * ~r 


2- 


iLi?  _ 4 M t n .2d 
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With  (his  ri'sult,  (11.22)  is  rewritten: 


w-i 

$>lxo,ijoJN,H')  = 4 t CwkCN  H)  t 2q  X jlN-D*  H)  • ^it^j 

+ Z - [l^"atl-".-‘  -y  4Vl]TW^,g.)  ipfc^'»WW.(«,|j.1»Hl) 
r (ij.z-0 

where  DWTE  is  the  sidewash  due  to  the1  wake  trailing  edge  mode  described 
in  section  iU.A.4),  with  spanwise  variation  n;  and  DWT  (x»,tja  ^N- H)  is  the 
sidewash  due  to  a triangular  distribution  of  vorticitv  (figure  II.  5)  having 
its  vertex  located  a distance  n-h  from  the  trailing  edge  of  the  sail  with 
spanwise  variation  n. 

Moving  all  known  quantities  to  the  right  hand  side  of  the  equation: 


r«, 


•) 


M-i 

= pl*,u.,N-H)  * 2 STXOTk.U1Jj.1(N-IkH)  Vt«KT 

* Z ( -.?—  * 1 ~^^S«kvlXoj^o)  1"  ^ f |Vf*  H ) 

(n.ix) 


2. 


Of  the  chordwise  modes  chosen,  only  the  first  mode  chordwise  makes 
any  contribution  to 


Therefore,  the  final  form  of  (II.  6)  is: 


Zi  ^ Cuk(W'H) •^WSh.h.lXojijo)  - •^2iCjijN,Hl*t>u)se,d>q|jd 

IYI  . . 

= PKj (N  Hl-  ZZ,  t>uT. H ) • 4kt  - 1 1 £&**, t K(t ^uj_ OqflJ 
-^^bure^u.^N-K)  in.n) 
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1) 

2) 


In  this  form,  the  solution  procedure  is  reduced  to  the  following: 
Choose  a set  of  control  points. 

Construct  (by  numerical  integration)  the  left  hand-side  of  (11.27) 
in  matrix  form . 


Construct  (by  numerical  integration)  a table  of  tt  ) 

and  for  each  spanwise  mode  at  each  control  point 

for  a sequence  of  steps  downstream. 


6) 

7) 


Construct  the  right-hand  side  of  (11.27)  from  previously  calculated 
values  of  f the  results  of  step  3 , and  numerical  integrations 

to  determine  ins  o*-  (*»r  at  the  cont  rol  points. 

Solve  the  resulting  (over-determined)  system  of  equations  by 
the  method  of  least  squares  for  the  unknown  values  of  Cim.*.  H ) 
Calculate  the  new  values  of  by  equation  (II . 23) . 

Increase  N by  1 and  return  to  step  4. 


Note  that  for  the  first  time  step, 
and  the  equation  (11.27)  takes  the  form: 


. becomes  an  unknown, 


V\AW. 


Hi 


6)  Force  and  Moment  Calculations 
The  sectional  force  and  moment  produced  by  each  chordwise 
mode  is  obtained  by  integrating  (II.  8)  across  the  chord  of  the  sail  section. 
Identifying  the  "unsteady  part  of  Ap  as  *he  ^ term 

and  the  "steady  part  as  the  term  one  obtains: 


■ - aS • ■J.-X.-, 


„ V - 


where:  S.  F.  denotes  the  unsteady  sectional  force  for  the  same  mode. 


Note  that  the  total  sectional  force  for  each  models  given  by  the  local 
sectional  mode  strength  times  the  steady  term  plus  the  rate  of  change  of 
the  local  sectional  mode  strength  with  distance  traveled  times  the  unsteady 
term.  These  sectional  forces  and  moments  are  then  integrated  (numerically) 
across  the  span  of  the  sail  to  yield  the  total  side  force,  yaw  moment,  and 
roll  moment  on  the  lifting  surface. 


« v «r.  \r*,r;'w  rvyw^ ; 
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The  sail  tip  streamline  trajectory  is  calculated  in  a step  by 
step  manner  from  the  following  relationship: 

d$(xl  htf (jj.j) 

ItvC-XeH^+^oViV^ 

^cu» 


r0n~  f 

n - 

X»*.  - Xqk-i  2w 


which  states  that  the  local  streamline  slope  is  approximately  equal  to 
the  local  radial  velocity  produced  by  the  source  distribution  representing 
the  hull,  divided  by  the  free  stream  velocity. 


B.  The  Interaction  Response  of  the  Hull  to  a Step  Change  in 

Sail  Ci  rculation 

The  simplest  means  of  analyzing  the  unsteady  fo  rce  s on  the  hull 
caused  by  the  velocity  field  of  the  sail  wake  following  a sudden  change  in 
sideslip  angle  is  to  consider  th°  response  of  the  hull  to  a step  change  in 
sail  circulation,  then  use  the  superposition  thee  rem  (convolution  integral). 

In  essence,  the  wake  produced  during  the  indicial  response  of  the  sail  is 
considered  to  be  constructed  from  the  wakes  due  to  an  infinity  of  infinitesimal 
step  changes  in  circulation. 

The  distribution  of  vort’eity  in  the  sail  wake  following  a step  change 
in  sail  circulation  is  as  shown  below: 

j ^ ^ Tip 


The  wake  resembles  a lifting  line  mo\ing  along  the  wake  trajectory, 
its  trailing  vortices  lengthening  at  a rate  equal  to  the  free  stream 
velocity.  For  simplicity  of  calculation,  the  trailing  vortices  were 
considered  to  be  piecewise  straight  as  shown  in  Figure  III.  4.  This 
simplification  should  introduce  little  error  since  submarine  hulls  ai“ 
typically  very  slender  bodies  (having  uniformly  small  slope  and  curvature). 


Figure  1111.4) 

Applying  '.he  law  of  Biot-Savart  to  a segment  of  the  trailing  vortex  sheet 
one  obtains  for  the  side  wash  velocity  at  the  hull  centerline: 


IS 


± 

ITT 


Where  z^  is  the  norma!  distance  from  the  trailing  vortex  sheet 
(which  may  lie  in  a plane  other  than  z 0)  to  the  point  Xq,  y^. 


u&Vi, 


Applying  Biot -Sava  rt  to  the  bound  element: 


Ar  < i ^ (*k-K>)  Ah.  _x 

*U"‘  = '^\  ' 

The  total  z component  of  the  velocity  due  to  the  wake  at  the 

instant  that  the  bound  element  has  travelled  to  x x is: 

n 

~\T  . - _ i f^rii  * tXfc-Xo)  Alv 
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The  other  two  components  of  velocity  due  to  the  wake,  V 


x wake' 


V , were  considered  to  be  negligably  small  compared  to  the  V 
y wake  K z total  wake 

component. 


•ait* : ■ i i ci-  i 


( . Force  and  Moment  Calculations 


The  hull  of  a yawed  submarine  is  typically  modeled  as  a distribu- 
tion of  sources  and  sinks  along  the  axis  of  the  body  to  produce  the  body's 
axial  shape  and  a distribution  of  dipoles  having  their  axes  pointed  in  a 


direciion  opposing  the  crossflow  to  satisfy  the  crossflow  tangency  boundary 
condition.  Both  of  these  distributions  can  experience  forces  and  moments 


due  to  the  external  flow. 

Following  reference  o)  . one  reaches  the  conclusion  that  the 

source- sink  distribution  will  experience  both  a force  and  moment  due  to 

the  presence  of  the  wake.  There  will  be  no  force  on  the  dipole  distribution 

sine  e the  sail  image  provides  for  no  influence  of  the  wake  velocity  field 

on  the  local  dipole  strength  required  to  satisfy  the  crossflow  boundary 

condition.  There  is  no  moment  on  the  dipole  distribution  due  to  the  wake 

was  a result  of  the  assumption  that  V , and  V , are  negligably 

y wake  x wake 


small  compared  to  V 


z.  wake 


The  force  on  the  source  distribution  as  a function  of  distance  travelled 
is  calculated  from  a relationship  given  by  McCreight  which  is  based  on 

slender  body  theory  and  Lagally's  theorem. 

1Sk.cw- 

$'U)uiX|S^<lx  ini.s} 


Where:  Wfx,  s)  is  the  local  sidewash  velocity  due  to  the  wake 


and  S'fx)  is  the  derivative  of  the  cross  sectional  area 


curve  for  the  hull. 
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Similarly,  the  moment  on  the  source  distribution  is  given  by: 

S lin. t) 


IV.  Rudder-Wake  Interaction 

A.  Following  the  philosophy  described  in  III.  B,  the  response  of  the 

rudder  to  the  velocity  field  of  the  sail  wake  will  be  obtained  by  finding 
the  response  of  the  rudder  to  the  wake  resulting  from  a step  change  in 
sail  circulation,  then  applying  the  supe  rposition  theorem. 

The  rudder  analysis  will  be  carried  out  on  a pseudo- steady  basis; 
i.  e.  , at  each  instant  of  time  the  rudder  an  its  wake  are  assumed  to  be  in 
a steady  state  condition  of  vorticity  distribution.  This  is  a reasonable 
assumption.  Typically,  a submarine  rudder  has  an  aspect  ratio  of  roughly 
the  same  magnitude  as  the  sail,  with  a chord  measuring  about  one-third 
of  that  of  the  sail.  Hence,  the  non-dimensional  dynamics  (transient  response 
of  total  bound  circulation  versus  number  of  chord  lengths  travelled)  of  the 
sail  and  the  rudder  will  be  very  nearly  equal,  but  the  rudder  will  respond 
roughly  three  times  as  quickly  in  real  time  due  to  its  shorter  chord 
length.  The  pseudo- steady  analysis  essentially  assumes  that  the  rudder 
responds  instantaneously  to  a change  in  angle  of  attack.  Since  tne  rudder 
responds  three  times  as  fast  as  the  sail,  the  pseudo-steady  analysis  should 
give  a fair  representation  of  the  interaction  of  the  sail  wake  and  the  rudder. 


B.  Rudder  Response  to  Step  Change  in  Sail  Circulation 


The  details  of  all  of  the  steps  required  to  analyze  the  response 


of  the  rudder  have  been  covered  in  previous  sections.  The  appropriate 


integral  equation  is  given  by  II.  10.  The  equivalent  span  is  calculated 


from  II.  18.  "'he  spamvise  and  chordwise  loading  modes  are  those  used 


for  the  sail  (the  special  chordwise  mode  is  not  used).  Equation  II.  10 


is  solved  by  choosirg  sixteen  control  points,  integrating  the  loading  modes 


numerically,  calculating  the  left  hand  side  (local  angle  of  attack  at  control 


point)  from  equation  III.  4,  and  solving  the  resulting  system  of  equations 


by  lea s t- squa  res . The  forces  and  moments  on  the  rudder  are  calculated 


in  the  same  manner  as  those  of  1 1 sail,  with  the  major  exception  that  only 


the  steadv  contributions  of  each  mode  are  considered. 


Supe  rpos  itio.  Theorem 


Consider  a two-dimensional  section  of  the  sail-wake  vorticity 


distribution  as  shown  in  Figure  V.  1. 


taU-T.j) 

ri  — 


R<yj*fc"V.  1 


The  wake  can  be  considered  to  be  constructed  from  an  infinity  of  wakos 


from  step  changes  in  sail  circulation  of  strength 


?*■-**--  , * - ^ I 
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From  relationship  II.  4: 


!SHx,u,->WlJt,eU-?.^  * - ^ (v.i) 

If,  for  example,  the  velocity  field,  and,  hence,  the  force  on  the  hull 
due  to  the  wake  of  a step  change  in  circulation  is  known,  the  force  on  the 
hull  due  to  the  wake  of  the  indicial  response  of  the  sail  can  be  calculated 
from  a convolution  integral'. 


= - ( r^jy)  ^ h otf  a? 


Where: 


) 

Y,  ,,  (s)  The  side  force  response  of  the  hull 
hull 

to  the  wake  produced  by  the  indicial 
response  of  the  sail 


^ull  ,3) 


The  side  force  response  of  the  hull 
to  the  wake  produced  by  a step  change 
in  sail  circulation 


r fx-T,  0)  The  instantaneous  mode  strength  of 

the  vorticity  being  shed  into  the  wake 


I ?.d  typically  has  very  large  derivatives  at  s 0,  so  that  from  a numerical 
standpoint  it  is  advantageous  to  integrate  V.  2 by  parts 


V«-  -Y^C?)rtttvT)  j’  . !-3I4t)I’„.Is-tUT 
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Similarly,  foi  (hr-  remaining  forces  ant!  moments  on  the  hull  and  sail: 


= N£jo)r**is) 
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VI.  Unsteady  Response  Results  for  Example  Vehicle 


A.  Vehicle  Description 


The  analyses  described  in  the  previous  sections  were  applied 
to  an  example  hull- sail  - rudde  r combination,  the  dimensions  of  which  are 
shown  in  Figure  VI.  1.  For  simplicity  of  calculation,  the  sail  and  rudder 

planforms  were  both  chosen  to  be  rectangular. 

— td  1<J—  1. 


The  bow  and  stem  portions  of  the  hull  were  assumed  to  be  parabolic 
i r«»  ) while  the  remainder  of  the  hull  was  approximated  by  a 

second  order  interpolating  polynomial  (ff  +CXV  ). 

B.  Indicial  Response  of  the  Sail 

The  strengths  of  the  vorticity  cross  modes  used  in  the  solution 
of  the  equation  (II.  M versus  the  number  of  chord  lengths  travelled  after 
the  step  change  in  sideslip  angle  are  plotted  in  Figure  VI.  2.  For  this 
aspect  ratio,  the  circulation  buildup  is  very  rapid;  after  the  sail  has 
travelled  only  one  chord  length,  the  lifting  cross  mode  C,^  has  risen  to 
roughly  85%  of  its  steady  state  value.  A rectangular  wing  of  aspect  ratio 
6 would  have  to  travel  roughly  4 half  chc ~ds  to  attain  a similar  percentage 


of  its  steady  state  circulation. 


Tlit'  .orces  and  moments  corresponding  to  the  circulation 


responses  of  Figure  Vi.  2 are  given  in  Figure  VI.  1 . a,  b and  c.  The 


force  and  moment  calculations  indicate  that  at  the  starting  instant  the 


sectional  forces  on  the  sail  are  increasing,  producing  the  initia1  "humps' 


in  the  responses.  The  presence  of  the  "humps"  is  a significant  result 


since,  for  example,  the  sail  will  operience  a maximum  roll  moment 


overshoot  of  28%  during  the  transient  response  of  the  sail. 


C.  Response  of  the  Hull  to  a Step  Change  in  Sideslip  Angle 


tSail-Hull  Interaction  Only) 


The  side  force  and  yaw  moment  on  the  hull  due  to  a step  change 


in  sideslip  angle  are  plotted  in  Figures  VI.  4.  a and  b.  Near  the  start 


of  the  response,  both  the  side  force  and  yaw  moment  are  positive, 


whereas  their  corresponding  steady  state  values  are  negative.  This 


effect  is  due  to  extensive  changes  in  the  velocity  field  due  to  the  wake 


that  occur  as  the  wake  lengthens.  Near  the  start  of  the  response  the 


wake  has  a velocity  (V  ) distribution  along  the  axis  of  the  hull  as  shown 


in  Figure  VI.  5.  a);  the  steady  state  distribution  of  V^  is  shown  in  Figure 


VI.  5.  b. 
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The  i ni ti a 1 distribution  of  wake  vorticity  produces  a negative 
on  the  forward  part  of  the  hull  (source  distribution)  and  a positive 
V on  the  after  part  of  the  hull  (sink  distribution)  and  hence,  bv  I.agally's 
theorem,  produces  a positive  side  force.  The  steady  wake's  predominant 
effect  is  to  produce  a negat’ve  velocity  on  the  after  part  of  the  hull  and, 
by  Tagally's  theorem,  a negative  side  force. 


D.  Response  of  Rudder  to  Step  Change  in  Sideslip  Angle 


The  side  force,  yaw  moment,  and  roll  moment  on  the  rudder  due 
to  the  sail  wake  are  shown  in  Figures  VII.  6.  a,  b and  e).  The  wake  has 
a comparatively  small  influence  on  the  rudder  until  the  trailing  edge  of 
the  wake  passes  by  the  rudder.  As  the  trailing  edge  of  the  wake  passes  the 
rudder  the  forces  and  moments  change  sign  (indicating  a change  in  sign 
of  the  sidewash  produced  by  the  wake  at  the  rudder)  and  rise  rapidly 
toward  their  final  values. 


E.  Total  Configuration  Response 


In  order  to  compare  the  relative  magnitudes  of  the  forces  developed 
on  the  sril,  hull,  and  rudder  and  the  relative  time  scales  involved  in 
their  responses  Figures  VI.  7.  a,  b and  c were  constructed.  These  plots 
represent  the  sums  of  the  total  transient  forces  and  moments  on  the  sail 
but  only  the  interaction  forces  and  moments  due  to  the  sail  wake  on  the 
hull  a i 1 rudder.  Sine'  the  sail  has  the  most  dominant  contribution  to  the 
force  and  moment  on  the  total  configuration,  reference  lines  for  the  steady 
state  forces  and  moments  or  :he  sail  were  added  to  Figures  VI.  7.  a.  barn 
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to  facilitate  a qua  nti  la  ti  vc  romparison  of  thr  importance  of  each 
interaction.  The  yaw  moment  was  referred  to  the  mid-chord  of  the 
sail;  the  roll  moment  was  taken  about  the  axis  uf  symmetry  of  the  hull. 

In  all  cases  iside  forc  e,  yaw  moment,  roll  moment)  the  initial 
portion  of  the  response  is  dominated  by  the  sail  response.  The  "hump" 
that  was  mentioned  in  Sec  tion  VI.  R indicates  that  initially  the  submarine 
acts  as  a lightly  damper]  force  and  moment  generator;  i.  e.  , a step  input 
of  sideslip  angle  causes  a rapid  rise-  high  overshoot  response  of  force 
anrl  moment. 

The  portion  of  the  response  between  one  and  five  chord  lengths 
travelled  is  charac  terized  by  the  interaction  response  of  the  hull. 

The  side  force  response  amounts  to  onlv  a small  percentage  of  the  steady 
force  or.  the  sail.  The  interaction  yaw  moment  of  the  hull,  however, 
causes  an  ultimate  moment  variation  of  20%  of  the  steady  moment  of  the 
sail. 

The  portion  of  the  total  configuration  response  past  five  chord 
lengths  is  characterized  by  an  abrupt  change  in  force  and  moment  due  to 
the  rudder-wake  interaction.  The  ultimate  change  in  side  force  caused 
by  the  rudder-wake  interaction  is  roughly  20%  of  the  steady  force  on  the 
sail.  The  yaw  moment  response  is  the  most  dramatic,  the  change  caused 
by  the  rudder  being  roughly  120%  of  the  steady  moment  on  the  sail.  This 
is  a result,  as  one  might  expect,  of  the  large  moment  arm  between  the 


£ 


mid-chord  of  the  sail  and  the  center-of-pressu  re  of  the  rudder.  The  roll 


Figure  VI.  3b 


Distance  Traveled  - Half  chords 


Side  Force  Response  of  Hull  to  . lradian 


dder  Interact 


Moment 


Figure  VI.  7a 

Side  Force  Response  of  Submarine  Hull -Sail- R udder  Combination  to  . lrad 


sos 


Distance  T raveled  - Halfchord  s of  Sail 
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Cnru’lus  1 on  s 

1)  Tho  transient  response  of  the  sail  to  a step  change  in  sideslip 
angle1  shows  a large  initial  overshoot  of  sideforce,  yaw  moment, 
and  roll  moment  with  re'spe-et  to  the  steady  state  values. 

2)  The  interaction  response'  of  the  hull  to  the  de-ve-loping  sail 
uake  shows  the-  development  of  a modest  side- force  and  a more 
size-able  stabilizing  yaw  moment. 

31  The  interaction  response  of  the  rudder  to  the  developing  sail  wake 
is  cha  ra ct  e*r i zed  by  a sudden  reversal  of  rudder  side  force  as  the 
trailing  edge  of  the'  sail  wake  passes  by  the  rudder.  This  force 
reversal  yeilds  a small  change  in  roll  moment  and  a very  large 
change  in  yaw  moment. 


>ftipyb-VL?  J^*.*  >-^  JP?-.^r'.y  ys^'j'.-y  ^ j-^i1  ^vj  - ♦y’^r^r- v ^ 4^  .^^-j- 


-«a& ^gai&»»t*^^ — ~ 1 .— -V..*  . . ' 1 ' ' ??.. 2 ?*  " 
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APPENDIX  A 

DERIVATION  OF  EQUATION  II.  1 


For  a distribution  of  dipoles: 


^1  0 i 


idLw 


Since  this  flow  satisfies  I^place's  equation,  the  integral  can  be  written: 


4- ft  TU,8'IL4^)*L(7)S  dx<^ 

Sai  I V 


CJcJm. 


Considering  a rectangular  wing,  the  boundary  conditions  on  t 
are  given  below: 
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Integrating  the  first  term  by  parts  in  the  x direction: 
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Integrating  by  parts  in  the  y direction- 
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A similar  integration  by  parts  proc  edure  yields: 
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Hence: 
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SAIL  STEP  RESPONSE  PROGRAM 


INPUT:  1)  Vehicle  Geometry 


2)  Side  slip  Angle 


3)  Trailing  Vortex  Sheet  Angle 


OUTPUT:  1)  Step  by  Step  Circulation  on  Sail 

2)  Step  by  Step  Forces  and  Moments  on  Sail 

3)  Exponential  Approximations  to  Circulation 
R esponse  of  Sail 


100 


POO 


20 


PI 


1 5 


30 


r; 

i.i 


G< 


l "V1  L 


PI 


1 * r 


OAT^  = 7M« 


PI / 5 5 /5  3 


P( >U ,? ) , “ws ( P0, l 3 i 


l < 2 C , 1 3 ) , P ( 20  ) , R 1 ( 1 2 ‘ . CUU(  5 , 40 ) 

1L  ( 20) ,WWSt2C  ,3)  ,CC( 12) ,TAU( 12 ,2) 


r ; mi  . c [ ( 

oim^ns  i r wj  (3,4  0) , wr, p(  3,20,40) 

0 IREN'S  I 0M  A (1  2 ) 

01  MENS  ION  USe  ( 12, 2U>  ,USM(  12 ,20)  ,'JSR  < 12 ,20)  ,TSF  (20)  ,TYM<20)  ,TRM(2C) 
0 I *■••■  Nc.  ! r-j  c r ( PO  , 1 3 ) , S r ( 4 ) , = M ( 3 ) 

P P A 0 1 U J , ' 4 vc  , N v $ , \'  C P C , N C P S » N S " ,H,AL,TP,0HS»FSl 
k ro  = 'jf  o<-  -*n(-  pc 


*v 

1/ 


nsps=nsp 

FORMAT  < SI  5 , 5PN.3  ) 

R2AD200, ( CP( I , 1 ) , 1 = 1,  NCP) 

FPCVAT  ( lS^  ) 

•N  •■•0  S = M ‘-  S ■*  \w  2 

CAICI'LAT-  > 0 U I V A L r-  Nj  T SPAN  C JUNCTION  CCCRDI  PATE 
AP=.5*(RSL-RHS)*(RSL+RHS)/PSL 
Y1  = AP-»  SL  +3HS 
YH=.6?66l«57*ARSIN(Yl/AR  ) 

CjTARL  ISH  S P A N ■>,  I Sr;  CCMPCL  PCTMS 
or  20  1 = 1 / roC 
C P ( 1 *-  ( 1-1  )*N0PR,2)=Yl  + .l*AR 

no  21  : =1 ,ncpc 
no  21  J=  7,n-CPS 

CP(Jf( 1-1 )*NCPS,2)=(2*J-3)*.5*AR/(NCPS-1) 

US-  2-'c(  - T H M,,  )-  7 0 SOLVr  PC  ON  FIRST  STEP 
HI = . 5 =H 
H2=l . HH1 
H3= 1.-H1 
AP1=AP/Hp 

no  1 1-1,  NCP 
XO=(CP(  ( , 1 ) -H 1 )/H? 

YO  = C 0 ( I ,.2  )/H? 

PW'l(I,N-irS  + l)=AL*(l.  + (RHS/(PSL-AP  + CP(l,2>)***2) 

DO  1 J=1,NMS 
YH1 =YH/H2 

WWS(I,J)=SINT(X0,Y0,J,0,AP1,YH1)/H2 
0V4S1  (1,1)  =WWS  ( I , J ) 
nn  is  1 =1  ,n<~  p 

X0  = CP(  1,1) 


YO  = CP ( I ,?  ) 


00  15  IJ=2,NMC 
K = ( I J-l  ) TNMS 
a0  15  J=l ,NVF 

OWS 1(1,K+J)=SINT(X0,Y0,J,IJ,AP,YH) 

N2=N3DS  +1 

PRINT50 

DO  30  1 = 1, NCP 

PRINTS  00,  (DWS1 ( I , J ) , J = 1 ,N2i 

TALL  GLSOt'TWS  l ,nl  ,1  L , NCP,  N"OS  ,RUC,,0  . ,0.  ) 

P 1 N ~ Tc  AIL  IMG  cCCr-  VAL'Jt  OF  VHRTICITY, 
PTPsSOk-  ( 1.  -(  h?/m?  )<■«?  ) 

V0»=H3/H2/H?/OI P 


ETC 


. M , 

!:# 


*?  ■ 


"TO 

t' 


■i 

1 


aui.vig.  .,  .n..  j,  ..... 


'•’AIN 


DATE 


74  IAS 


21/55/53 


I 


Sr 


I 


E 


I 


I 


I ! 


it 
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DO  2 J=1,N«S 

coin  J, 1 ) = -t .*81 (J )*DI P/3.1416 

2 CWJ( J ,1 ) =-l . *81 (J )*VOR 
NO  = NMS  +1 

no  I = NO  f Nwrs 
1 7 COU( I , 1) =81 II) 

F I IV  0 DOWNWASH  OF  UNIT  PSEUDO  MODE 
00  3 I =1  , NC  P 
XO=OP(  1,1) 

YO=C  P ( 1,2  ) 
no  3 J =1  , NM  S 

3 WWS(  I ,J)=(WWS(  I ,J)+(rH  P/3.1416)*(  SINT(XO,YO,J,  1,4R,  YH ) -Tf  OvC < A P , H , 
1 1 ,X0 ,Y0, J,TH,YH) ) J/VCR 

DO  4 I =1  , NC  P 

4 PRINT500, (WWS( l , J ) , J=  1 , NM  S> 

500  FORMAT ( J 0 * , 10X,10E12.4) 

00  5 K=l,N”S 
00  5 1=1, NCP 
XO=CP<  I , t ) 

1 — ro=CPii,zi 
DO  5 J=1,NSD 

5 WGF(  K , I , J)  = WAKE (H,AR, J,XO,YO,K,TH,YH) 

FIX  MATRIX  fouation  array  for  SUBSEQUENT  VALUFS  OF  1ST  HCDE 
DO  6 I =1 , NC  p 
XO=CP<  1,1) 

YO*CPt  1 , 2 T 
DO  6 J=1 , NMS 

6 DWS( I ,J ) = SI NT ( XO,YO ,J , 1 ,AR ,YH)— WWS( I ,J)*6.2832/H 
IFINMC.EQ.l IRQ  TO  13 

DO  12  1=1  , NCP 
XO=CP( 1,1) 

YC^CPtr,?!  — 

DO  12  J=2,NMC 
N=NMS 

DO  12  K=1,NMS 
L =N*< J-l) +K 

12  DWS  ( I »L)  = SINT (XO,YO,K,J ,AR,YH) 

r3-CONTTWUF  ~ 

DO  11  N=2,NSD 
DO  7 1 = 1, NC.' 

P(  D =0. 

XO=CP< 1,1) 

YO=CP ( 1,2 ) 

OCT'T"J^r,NMS  - * 

FIND  NEW  DOWNWASH  OF  TRAILING  EDGE  MODE 
BCI) =B  ( I ) +CCU ( J, 1 )*TECMC( AR,H,N, XO, YO, J,TH, YH ) 

F I NO  DOWNWASH  OIJE  TO  PL  MODES  IN  WAKF 
N1=N-1 

DO  fl  I J = l , N 1 

- n Rm^Etii+cnji j,i jfTvR^nrt — — 


■ . , .■  - i - - . i. 


LrVFl  21 


MAIN 


DATE  = 79  199 


5S< 

21/55/53 


7 8(1  )=P(  I ) - ( C W J ( J , N-  1 )+6.2832*C0U( J, N-l  ) /H ) * WwS  ! I , J ) 

DO  9 1=1, NCP 

DWSL ( I ,NMDS  + 1 )=B!  I )+AL*(l  .M  RHS/  (RSL-AR*CP<  I,  2 ) 3 1**2) 

DO  9 J =1 1 NMDS 

9 o9i  < i * j ) =rws  ( I * J ) 

CALL  GLSOIPWS 1 ,81 ,! L, NCP, NMDS  ?UG,0.,0.» 

PR  I NIT  5 0 

N3='NMUS  + 2 

DO  31  1=1, NCP 

DWS ( I , N3 !=DWS ( I * 1 1 *B1 ! 1 ) 

DWS  ( I ,N2)=P(I)+AL*(l.'»-(PHS/(RSl-AP4-CP(I»?)i)**2) 

DO  32  J=2,NM0S 

3 2 DWS ( I ,N3 )=DWS( I ,N3  >+Bl ( J !*DWS<  I ,J  ) 

"JI  TPTNT5CK) } LOWS! I,J1»J=1,N31 
DO  25  I = l »NMD S 
25  COU( I , N) =B1 ( I ) 

DO  11  J=l,NM? 

11  Cw J ( J , N)= (6 .2832/H) *< COU( J ,N )-COU( J , N- 1 ) ) -C W J ( J , N - 1 ) 

00  23  1=1, NCP 

' xo=cpri,n  ' 

YO=CP( 1,2 ) 

ST(  I , NMS  + 1 1 =AL*(1.MRHS/(RSL-AR+Y0)  )**2) 

DC  23  J=l,NMS 

23  ST(  I ,J)=SINT(XO,YO, J,0, AR,YH) 

CALL  GLS0(ST,B1 ,1 L,NCP,NMS,BUG,0.,0.) 

CAL.U  _STTTffRVTH,TlHS  , S F , R N ) 

PRINT  800,  AP. 

800  FORMAT!  »1*, 10X,*STFP  RESPONSE  OF  SUBMARINE  SAIL  WITH  HULL  IMEPFER 
1 ENCE  , A R=  ',F5.8) 

PR  I NT  400 , H 

900  FORMAT ( 'O' , 10X,  ' 8 E T A = .1  H=',F6.9,'  HALF  CHOPCS') 

P1HRT&75  

625  FORMAT! *0* , 10X , ‘STARTING  MODE  STRENGTHS'! 

1 1 = 0 

PRTNT6  00, (I  I , I , PI ( F) , 1 = 1 ,NMS  ) 

DO  10  N= 1 , N SO 
PR INT3J0.N 

~-30U  FORMAT!' O',  IOXF'WATCE  VENGTH=' ,13, 2X, 'STEPS'! 

DO  10  I = l , N MC 

10  PRINT6  00, ( I , J , COU  < < l-l)*NMS+J,N), J=1,NMS1 
6 OC  FORMAT! IX  ,3! ?X, *8! ' ,1 1 , I)  ,'  )=' ,6 12. A) ) 

PR  I N T°0O 

900  FORMAT  ( • 0 • ,10X,'WAKT  MCOf  STRENGTHS') 

DO  19  J=1,MMS 

19  PRINT  700,  (CWJ(J,N},N=l,NiDI 
700  FORMAT < '0e , 10X.10E12. 9 > 

'•.ALL  STrAOY(AR,AL,TH,NMC,NCP,rP,YH,PHS,PSL,CO,CWS,NMSl 
P R 1 N T 8 5 0 

850  FORMAT  ('  U*  ,10X  STL  l n\  STATE  Mrpr  STCFNGTFS') 

00  18  1 = 1,  NMC  •-*- 
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L c V L "1 


!f  T : = 7 4 14  9 


21/55/53 


l'1  rt  i,  -/  ,(!,!,  ( ( i-l  )~T  R + J)  ,J=l,'  yS) 

r PALI  7 A"R  <K  vr  ,MV«5  , MSr  , CC'J  . CO,  h,T  AU  . A ) 

p c i • i j 5 s o 

559  R n*-'  *.t  ( i ■)  ■ , iox  . • R XP°\tNT  l AL  APPROXIMATIONS  TO  *-005  HTSTCRIES*  ) 
f 1 . !=1-’  ■■■r 
' ( If  1 = 1 , 
k = j-c  *(1-1)0 
A 1 1 .-A! < ) 

L 6 PF  I NT  »’  50,  I , J ,00  < K ) , A l , TAU(K,  1)  , A ( K)  , TSUI  K ,2  ) 

750  f HRMAT ( • 0*  , IOX  , '3  ( • , 1 1 , •-• , 1 1 • ) = ' , r 12.4,  ' ( 1.-  • , F 1 2 . 4 , * c XP ( • , E 12. 
It,  ' * S ) - •,-12.4,'RX‘>(',::12.4,'  = S))') 

CALL  ^S,'\Fr,cSLtivHS,AP,CrJU,5C,:WJ,H,USr,USM,:j<:P,T?FfTVP,TP 

1 M) 

T”  PR7NT910 

“■  910  FORMAT ( 'O' , IOX,  'FORTES  5 MOMENTS  OF  MOOES') 

R PRINT920, (USF ( 1 , J)  , J = 1 ,NSOS ) 

PPIN'r930 

920  FORMAT  { «9*  , IOX,  10212. A ) 

930  FORMAT ( ’O’) 

EJ  ,7"  ,-R  1 NT920 » (USM(1  , J ) , J=  1 , NSDS  ) 

rcl  pr  i nt  930 

‘■'•j-  PRINT92G,  (USR  ( 1 , J i , J=1  ,NSOS) 

^0  2 2 1 = 1 ,’JMS 
PP  1NT930 

2 2 PR  INT920,  ( USF ( I+NMS , J ) ,J  = 1,NS3S) 

j PR  j NT940  

\ 940  FORMAV(  'O'  , IOX, 'TOTAL  SIDE  FORCE') 

J.  PPIM1920,  (T  S F ( J ) , J=  l , NSDS  ) 

PRINT  950 

950  FORMAT!  'O'  , IOX,  'TOTAL  YAW  MC Mc  N T ' ) 

! !•  PP1NT920,  ITYM(J),  J=1,NSDS  ) 

tar PRINT  960 

960  FORMAT! *0* , IOX, 'TOTAL  ROLL  MOMENT') 

F.  PRINT920,  (TRM(J),  J = 1,NS0S) 

T CALCULATE  STARTING  FORCES  C MOMENTS 

SSF  = - 1 .5  70R«R1 ( 1 ) *SR ( 1 ) *2. 

SYM=  SS  F*  . 5 

fP ':~^SRN=-RFtT  ll*SSF/SFt  1 ) * 

s DO  24  1=2, NMS 

l SSFA=-1 .5708*B1  ( I )*SF  ( I 1*2. 

1 ‘ " SSF=  SSF  + S S R A 

SYM=SYM+. 5*SSFA 
24  spm  = SPm+dm(  T ) *SSFA/SF  ( I ) 

F -SRM=SFM^SF  *(RSL-AP) 

PP  INT650 

650  FORMAT!  'O' , IOX, 'START ING  FORCES  C MOMEMTS') 

11  PRINTF.75,  SSR,  SYM,  SRM 

675  FORMAT | • O'  , 2 X , ' S F = • ,F12.4,2X,'YM=',f12.4,2X,'PM=',F12.4) 

C CALCULATE  STEADY  FORCES  £ iC'FMS 

ESF=3.1416*CD! 1 )*SE!1 )*2,  ' " ~ 
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*~‘~  'Jf&tm&JaHe H. 1 _ 
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01  tC 

TITBIT 

0189 

0190 

■Ji9  r 


oir-r 

019  3 

ttpit 

0195 

0196 
0T  97 
0198 


0199 

I uzao  - 

0201 

0202 


V 


9 Y " = . 5 * ? S r 

FFM=ESr*RM<  1) /SF { 11 
on  26  1=2, NMS 

ESFA=3.1416*CC( I)*SF( I )*2. 

e sp=r  s f*  sp  a 

F Y ->1=  f Y *'  + . 5 * : 5 F A 

26  f 9 M = lR  M+  SF  A*  P^‘  ( I ) / S 9 ( I ) 

EPM  = FFM+ESF*(  RSL-AR) 

DO  27  1=1 , NMS 

27  EYM=EYM+. 7854*2. *SF (I )*CQ(NMS 
PR INT  660 

660  FORMAT  CO*  . 10X, ‘STcAPY  STATE 
PPINT675,r‘SF,EYM,F2w 
50  FORMAT  Cl*) 

STOP 

END 
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r 
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-pel  t i\m  IV  0 L : V P L ? 1 


SINT 


j 

0 00  1 

rirrur.N  r!NT(x,Y,N,u,A=>,Y‘  ) 

I w** 

— DW2 

DIMENSION  Rf6),F(fc),W(3),F(6) 

1 

0003 

Z 1 ( A , B ) = B* ( l.  + A 1-1. 

I 

i . 

0004 

Z 2 ( A , P ) = R* ( 1 . -A  )+  A 

“otms 

Z 3 ( A , P ) - 1 * ( 1 . +4  )+  3 

1 

000  6 

Z4  ( A , P ) = 3 * ( 1 . -A ) ♦ 1 . 

0007 

ZTTA  = . 63  jhl  )7*A  Rc  If  ( Y/  A8) 

— — — 

007J8 

R ( 1 ) = .2386192 

0009 

R(2)=-l.*R<ll 

0010 

R 1 3) =. 6612094 

00 1 1 

P(4)=-1.*R(3) 

0 01  2 

p ( 5)  = . 73  246  8 6 

001  i 

P ( 6)  =- 1 . *R  ( 6) 

f . 

~ 0TJT4-  ~ 

W(l)= .4679139 

001  5 

W (2) =.3607616 

0016 

W < 3)  = . 171  3245 

7)017 

IF(7RTA.LT.O. ) Gr  To  3 

0018 

0 1 = 7 ETA 

0019 

02=1 .-ZFT A 

0027) 

00  1 1=1,6 

002  1 

A 1 =P  ( l 1 

L - 

0022 

E ( I ) = Z 1 C A 1 , ZET A ) 

0023 

F (I) =Z2( A 1 , ZETA  1 

0024 

P 1 = F ( 1 ) 

0 02  S 

P 2 = F ( 1 1 

i 

00  Z6 

FTAl=AR*S  INC  1.5708*61  > -Y 

i 

0027 

ETA2=A°*S IN ( 1 .5708*82 )-Y 

002  8 

Cl ( I ) =CINT( X.ETAl  ,M) 

0 02  7 

1 i?( I )=CINT{ X,cTA2,m) 

00^  0 

Gf  1 TO  8 

0 0 ' 1 

3 01=1 .+ZETA 

r~ 

i 

003  2 

D2=- I .*ZETA 

0033 

00  4 1=1,6 

i 

0034 

A 1 =R  ( I ) 

— — 

0035 

F ( I ) = Z 3 ( A 1 , Z E T A ) 

0 03  6 

F ( I ) =Z4<  A1 , ZE  TA  ) 

003  7 

P l = E ( I ) 

r 

CTU38 

R2=F ( I l 

! 

003  9 

ETA1 =AR*S  INC  1 .5  708*B1  ) -Y 

0040 

ETA2  = AR*S INC  1.570 8*0 2 )-Y 

004T 

Cl  (I  ) =CINT(X,ETA1 ,M) 

) 0 4 2 

4 c 2 ( I ) = C I N T ( X , F T A 2 , M 1 

0)4  3 

5 01-0. 

I 

~UTJ4TV 

G2  =0. 

1 

004  8 

on  2 11=1,3 

• 

0046 

•JJ  = 2*  I I 

L. 

0047 

1 J = 2 J-  I 1-1 

0046 

X 1 = F ( I J ) 

0 04  9 

X2  = r u i) 

r"~ 

OU50 

Y 1 =F  C I J) 

'IT! 


,C1(6)  ,C2(6) 


7 4 K9 


r 


/£ 


r^  1. 1 Wfr.  tY^V.’Ti 


es< 


nHH4N  IV  0 ICV^L  21 


OATF  = 7 '<14-, 


0 05  1 

0052 

0053 

0054 
005  5 
0 05  6 
005  7 
005  BT 
005  9 
0060 


Y?=f  ( J J) 

W1  = C1<  I J ) 

W2  = C 1 ( JJ J 
P l=f  2(  I J I 
R2=C2( JO ) 

01  =01  *W  ( m *(  n (XI  f Y,W1  , AR  ,N,Yh)  + F1  ( X2,  Y,  W2.4K  ffj 

f-2=G2+W(  II)  *(  F 1(Y1 , Y,  R1  ,AR  fN,  YH)  + F1  { Y2,  Y,  R2  , A«  ,N 

S I NT  = ( 01  *D1 +02*02  ) *. 19635 

RETURN 

END 


,YH)| 
, Yh  ) ) 


1 jjgwr 


triMOs&fo  :iS&cn±jikiSi  ±!stfe  HatHiki  iiiiSiu: 


f y ’y**  ■sy ’vyi* *v.’ *-> Ny* ,Ly.i 1 fiy^' ». 


IV 


ilV  " I ) ■ I V ( X , ' T ;'.  , 1 ) 

01  virus  KIN  R(6),E(6),F(6),W(3) 
ZllA,FJ=4*(l.*B)tR 
Z 2 ( A ,B) =3*  ( 1.  -A )+  1 . 

ZM A,  F >=  3M  1 . +A )-l  . 

Z4<  A , F ) = JM  1 . -A  ) + A 
I =L=.6366197*4PSIN<X) 

R ( 1 ) = . 2336192 
0(2) =-l.*P(l ) 

° ( 3)=. 6612094 
R(4)=-1.*P(’] 

P ( ^ ) =.9324695 
R (M  = -1.=P  ( 5 i 
W(  1 )=.  46791  30 
W ( 2) =. 3607616 
W<  3 ) = . 171  3245 
I F ( PF  L . GT . 0 . )GP  TF  10 
on  i i=i,6 
Y =R  m 

E(  H =Z  KY.DELI 
F(  I)  = Z2(Y,0EL) 

A1=1.*0EL 
A 2=- l.*nsi 
on  tij  20 
no  2 I =1  ,6 

Y=R  ( 1 ) 

E ( I ) =Z  3 ( Y , 0 EL  ) 

F ( I ) =Z4(  Y ,DEL  ) 

a l =PC L 
A2=l .-PEL 
r INT=0. 

DO  3 1=1,3 
X 1 =E  ( 2+  I -1 ) 

X 2=E ( 2*  1 I 
Y1  = F(?*T-1  ) 

Y2=R( 2«I ) 

0=A1*(  FF(  XI ,FTA,X,M)+FC ( X2,ETA.X,M)  J 
C=A2*<FCfYl ,ETA,XfM)+FC(Y2,ETA,X  M)  ) 
C I NT  =C  I NT-*- W ( I )*  (P*C) 

RETURN 

Fwn 


OA  T f = ?41  /,9 
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t rV:L 


21 


t-r 


iA< 

nzr-.  - 7 M a -7 


O'M.  1 

0002 
00C3 
0004 
30u  5 
0006 
000  7 
0003 
0009 


r i j-Nir.  t t ' tnx,':TA,xii ,«) 

IFCM.F0.0>G=S!N( 1.5 708* XI 
IF(M.EQ.l  )0  = 1 .-S  I N ( 1.5708 -XX  i 

I F ( M • r 0 . 2)0,  = -.5l»(CCSC  3.  L416*X)tSIN<  1 .5709*X  ) > 
Ir(”'<.rT.2)O=.5*(CCr)(,-'*1.5700  + (l.+X))*COS<(1'-?) 
A=9I\'(  1 . 9 708*  X J 

FC  = G*SOPT((A-XO)’>(A-XO)+:TA*£TA)/(A-XC) 

RETURN 

EN° 


. 5',C8*  ( 1 . + X 


1 


r" 


i 

.j 


F 1 


)Mr 


7414R 


C5< 


P'  - T C ft'l  1 V 


l V l 21 


r 
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0 0 0 1 

FUMC.T  J "V  PI  (A  ,V  , r,  AP,  N,  YH  1 

0002 

IFf  N.FO. 1 IOC  Tn  l 

0003 

IF(N.EQ.2)GC  TO  2 

0004 

IF  ( N.  £ 0.  3 ) 00  TP  3 

000  5 

1 

CONTIN'F 

00  of- 

IF(/.[T.YW)0P=5  IP  ( 1 .5703*  ( 7 + 1 . ) / ( YH+  1.  / > / ( YH+  1 . ) 

000  7 

IP  </  .Or.  YH)  r,S>-r  rz  (1  .5  700*  ( 4-1  . ) / ( YF-  1 . ) ) / ( YH-  1 . 1 

0008 

00  TO  4 

oooo 

? 

PONT  I NUF 

0010 

1 F 1 A . LT.  YH)OS=r  OS  ( l . 570  8*  ( A*)  . 1 / € YH*  1-1  ) / C YH-t- 1 . ) 

301  1 

1 r ( A • GF . YH ) G5  = Fr  c ( 1 .C>7C8*A/Vh  ) / YH 

0012 

! p ( A . r,P  .0.  1 GS  = -?  .*POP(  A * 3 . 1 4 1 6 ) 

001  3 

G(  * i 4 

0014 

3 

CONTI  Mjr 

0015 

IF  ( A.LT.  YH)  GS=3.*C  OS  (4. 7124*  (A+1.)/(YF  + 1.J1/(YMI.) 

00  i 6 

lr« A .Gt.YH)GS  = 3.*C0S( 4. 7124*( A- 1. ) /( YH-1.  M /(  YF-1. ) 

001  7 

4 

n=P*r,S/(AD*SINU*l  .57O0J-Y) 

001  8 

Pf  TURN 

001<? 

FNO 

L.. 


.r 


f!(>< 

" l T t(  A f" 

IV  0 l-vc  L 

2 1 W A k 

OAT  F = 74  149 

O'lt  1 

FlJ’Jf  T } • ’ '-J  WAK  (H  ,A°  , XO,  YO, 

I I , TH , YH  1 

0002 

DIMENSION  R (6 1 , E (61 , F(6 1, W( 

31 

r?." 

.<yc 

0003 

Z1 ( A ,B ) = 3*( 1 . +A  )-l  . 

r/ 

0004 

Z2<  A ,B 1 =«*(  1 .-A  1 + A 

0005 

Z3  ( A , •* ) = **(1  .+«)tP 

0006 

Z4(  A ,8  ) =a*  ( l.-A  1 + l . 

0 00  7 

ZrTA  = .6366197*ARSIN(Y0/AR  ) 

ctoob 

R ( l) =.2386! 92 

■ *■  * 

0009 

R(2)=-l.*R(  n 

0010 

R(31=. 6612094 

001 1 

R ( 4 1 =-  1 . * R ( 3 1 

0012 

H ( 5 i= .9324695 

001  3 

P ( 6 ) = - 1 . * P.  ( 5 1 

~00T«r 

W ( 1)=. 46791 39 

*■* 

0015 

W( 2)= . 3607616 

0016 

W( 31 =. 1713245 

0017 

IFIZFTA.LT.O. 1G0  TO  4 

001  8 

Dl  = Z ET  t 

001  9 

02=1 .-ZFTA 

VOZV  - 



DO  2 1=1,6 

• -T~r 

hi 

0021 

A1=R(  I 1 

'* 

0022 

E( I) =Z l< Al,ZETA  » 

0023 

2 

F(I)=Z2iJ 1 * 2 2 T A ? 

0024 

GO  TO  6 

002  5 

4 

n i=i.*z=tA 

0026 

D2  = -I .*7  ETA 

. . 

002  7 

DO  5 I =1 ,6 

0028 

A 1 = R ( I ) 

0 02  9 

R (I)  =Z3<  41  , Z R T 4.  1 

0030 

5 

F ( I ) = Z 4 ( A 1 , ZETA  ) 

003  1 

6 

W A K E = 0 . 

-OUST 

DO  3 I =1  ,3 

— • 

0033 

R 1=  E ( 2*1-11 

0034 

R2=  E (2*1 1 

003  5 

FI =F (2*1-1  1 

0036 

F2  = F( 2*1  ) 

00*7 

G 1 = 0 1 * ( f 9 4 ( o 1 ,/p,h,XO,YO,N, 

I I,TH,  YH1+FWA  (R2,AP.  ,F,X0,YC,N 

, ii 

1 1 

'T  ■■ 

0038 

G2=02*  (F'9A(  FI  ,AR,  H,  XO,YO,N, 
1 ! 

I1,TH,YH)+FWA(E2,AP,H,X0,Y0,N 

til 

0039 

3 

WAKF  = WA‘<'=4-U(  I 1 * ( G 1 + G 2 1 

0 04  0 

RRTijP  N 

004  1 

f N 0 

€7< 

1 - - \!  ! V ' I V ■ L 7 l C*A  0&TC  = 7 ^ 1 ^ S 


r 1 lK,f  t I , F-h  A(  Y,At  ,H,XO,YO,N,  I I ,T,YH) 

GS  = PS  P AN  ( Y » I I * Y H I 
XNl  = --j*H-< XO-l . ) *CCS  (T  ) 

XN2=XN 1 +H 
XN3  = XM  -H 

Z 0 = A fl  S ( X 1 - 1 • ) * S I N ( T ) 

Yl=Au*SIM(1.570fi*Y)-Y0 
Y2=A3M  Y1  ) 

S = Y 1 / Y 2 
PAft=Y2/ZJ 

S»l=SC=’’r  (Xf'.'l^XNl  + Yl^Yl+ZO^ZO) 

SR  ?=  SC*T ( XN 2*  XN  ?*  Y 1*  Y l+ZO* 7 0 > 

SR^SQPT  ( XN3*XN?4Y1*Y  l-»ZO*ZO  ) 

SI =T  QP  T ( XN1 #XN1 + 70*  ZO ) 

S2=SQFT(XN2*XN2+Z0*Z0) 

S3  = SORT ( X N? *X N3 +ZO* ZO  ) 

A=Yl/( Yl* Y1 +Z0*  ZO) 

G1=A*{  .5*  (XN2*SF2+XN3*SR3  )-XNl*SP.  1) 

02  =Y1  ^ ( AL  rr,(  XM  +S«1  )-.5  *<ALnC-{XN2+Sa  2 H-ALOG  (XN3  + SR3  ))  ) 

G3=  S* ( XN2*ALOG(  ( S P 2-Y2  ) /S  2 ) + XN3* A LOG  ( (S  R3-Y2  ) / S3  )-2  .*XN1*AL 
1 -Y2 I / S 1 I ) 

G4=Z0*  S*(  ATAN(RAR*XN2/S«2  1+-ATAN  (RAB*XN3/SR3  )-2  .*AT  AN<  P-AB*XN 
PWA=(  G1+G2  + G3+G4) *GS*.  125/H 
FWA=FWA*CCS (T I 
F’CTURN 


Uv  1 

— 0002 

0003 

0004 
OOUS 
000*5 
000  7 
OOOB 

0009 

0010 

0 31  1 
0012 

001  3 
“OOT4 

001  5 
001  6 

0017 

0018 
001  Q 

UOZO 


0-021 

U077~ 
002  3 
0024 


- . , r-«jn  n.’r*a  tmtc+rW*  n\jn*^~-***xv** 
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p . ' 4 \i~  iv  r. 

L f:VpL 

21  TCrM[l 

out  I 

c U‘1C  T I IN  Tf CMC! AP ,HtU,X0,Y 

0002 

DIMENSION  R(6),WI3),F(6),F 

1 

0003 

Zi(A,B)=B*(l.+A)-l. 

> 

0004 

Z 2 ( A ,P) =B*<  l.-Al  + A 

0005 

Z 3(  A * P ) = A*<  l.+BH-P 

0 00  6 

Z4(A,P)  = 3*(  l.-A  )H  . 

0007 

ZETA =. 63661 S7*A RSI N ( Y0/ AR ) 

__ 

— 0008 

R ( 1 ) = . 2386192 

1 

0009 

P(2)=-l.*R(l) 

0010 

R ( 3)  = . 6612054 

001  1 

; (4>=-l.*R(3) 

001  2 

F ( 5) =.  9324655 

0013 

R(6»  = -1.*M  5) 

r 

00T4 

Will  =.4679139 

i 

0015 

W(  2)= . 3607616 

0016 

W(3)=. 1713245 

001  7 

1F(ZETA.LT. Oi )GC  TP  4 

0013 

D1=ZLTA 

001  9 

02=1 .-ZET  A 

i 

OTJ7TT 

DO  2 T =1 » 6 

i 

0021 

Al=R ( 1 1 

i 

0022 

E ( I ) =Z1 ( A1 , _FTA ) 

*“  ■ 

002  3 

2 

F(I)=Z2(A1»7ETA) 

0024 

GO  TO  6 

0025 

4 

m =1 . + zeta 

1 

- TTO 2*6 

’ 

DZ=-1 .*ZETA 

1 

0027 

00  5 1=1,6 

0028 

A 1 =P  ( I I 

0029 

F ( 1 1 — Z 3 ( Al, ZETA) 

0030 

5 

Fill =Z4( Al , ZETA  ) 

0031 

6 

T ECMO  = 0 . 

003? 

DO  3 1=1,3 

• 

0033 

Y 1 =E ( 2*1-1 ) 

0034 

Y2=F<  2*1 ) 



“ TT075 

FI =P 12*1-1 ) 

0036 

p2=F( 2*1 ) 

0037 

01=ni*(PT(Yl,AR,H,X0,Y0,N, 

— 

OOT8 

G 2 =D  2 * ( F T (PI fAR  »H,XO»YO,N» 

0039 

3 

T EC,MD  = TEr  MO  *W(  I1*(G1  + G2» 

0040 

RETURN 

0041 

END 

— 
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OATF  = 


cs< 

7 4 14  3 


■T*'~ 


"T&i 

5 


i4  ! 


i vc-1 


M 


i:s< 


’V  ••  l L V’  1 ?\ 


'TT  = 7 ^ 1 4Q 


2 1/85, 


OOul 

rpoc? 

000  3 
0004 

WO  5 
0006 
0007 


0010 
001 1 
001  2 
001’ 
-ours 

00(1.5 

0016 

U0T7 

0018 

0019 


F'!rT!  ■’■  FT  (Y  , AR  , H,  XJ,  YO,  N,  ! 1 , T , YH  ) 

GS=F5P  AM  V,  1 1 ,YF) 

XN 1 = N $H— ( XO— l.)*CQS(T) 

X N ? = X N 1 -H 
XN3=XM?-H 

Yl=AR*-$  INI  l .S70R*Y)-Y0 
Y2=A9S ( Y 1 ) 

S=YL/Y? 

Z 0 = A 9 S ( X 0 — 1 .)*SIMT) 

PA9=Y2/Z0 

$»  1 = SCjC  T(  XN  1*XN  U Yl«Yl  «-ZG*ZC) 

SP.?  = ROK  T ( XN2*XN2+Y1*Y  1+Z0*70  ) 

51  =SOF  T ( XN1*XM+Z0»  70) 

S2=  SQF  T ( XN2*XN2+Z  0*ZO) 

A = Y1/ < Y1*YI+Z0*Z0) 

Gl  = . 0625v  (1./<XN1*XN1  + ZC*Z0)  + 1./(Y1*Y1  + Z0*Z0)MXM*Y1/5R1 
G2=AMXN?±5R2-XN3*SR  1 ) 

G3=Y1*M_-Gi  (XM+Sfil  )/  (XN2+SR2  ) ) 

G4=?.«'XN?"S*AL0G(  Sl*(  SP2-Y2)  / ( S 2*  ( SR  1-Y?)  ) ) 

G5=  Z0*S*  ( AT AN(RAB*XN2/SP2 )-AT AN (R AP*  XN 1 / SP 1 ) ) * 2 . 

FT  = 1.57U8MGl  + (.l2  5/(  H*H)  M < G2  + G3 + GMC5  ) ) *G  S/  ( 1 . - . 5 *H  ) 

FT  = FT #C OS ( T ) 

RFTURN 
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21  8 Tf  A n Y OATF  = 

)J)1 

S!<<-  l,r|J  ^7'  £OY  ( ftj  , ,'L  tTF,  N,JC  , NCP,CP,  YHf  »HS, 

0002 

ni  Mf  NS  ION  C PC  20, 2 ) , CHI  17)  ,CWS  (20, 13  > ,IL(  20) 

000  3 

DIMENSION  5(20,13) 

0004 

NNOS  = NMC  *NM S 

000  5 

rn  1 I = 1 , Nr f 

0006 

X0  = ( ■>(!,!) 

00t.  7 

Y 0 - r <J  ( 1,2  ) 

0008 

0WS(  I ,N*ns+l  ) =A!  * ( 1 .4-  (RHS/  CRSL-AR+YO)  )**2  ) 

0000 

OCi  1 J = 1 , NM  S 

0010 

l 

DWS  ( I , J 1 = S I NT (X0,Y0»J,1»AR»YH)-TRL(AR,YH,J.X0 

001  1 

iN  1 =N Mf  S + l 

001  2 

N2=NNHSf? 

0013 

on  2 1=1, NC  0 

~~  0014 

on  2 J = 1 , N 1 

0015 

2 

s(i,j)=r^s(i,j) 

0016 

CALL  GLSO(DWS ,C0, I L ,NCP,NMDS ,EUG,0. ,0. ) 

0017 

PR  INT50 

001  a 

m ? I =1  , NC  P 

001  5 

S(  1 ,N2!=M  I ,1  )*C(K1) 

otrzo  — 

DO  4 J=2 , NMDS 

0021 

4 

S(  I * N 2 ) = S ( I , N2  ) +5  ( I ,.J)*COUI 

0022 

3 

PRINT  SOO, ( S ( I , J ) » J = 1 ,N2 ) 

0023 

80 

FORMAT  ( » 1 * ) 

0024 

500 
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0025 
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“ 0 00  9 
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0007 
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001  0 
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001  J, 
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r 001  9 
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y ooi6 

W 3) =. 1713245 

^ 001  7 

I o (r  -l  .r,T  . j . ) y.  ti 

001  4 

R7  1 1=1,6 

0019 

Y = c ( ! ) 

nrwzu 

P(  I)  =Zl (Y,9FL ) 

f,  0021 
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p U ) =Z2<  Y.OFL ) 

| 002  2 
0023 

A l = 1 ,+m-l 
A 2 = - 1 . *0  = 1 

0 0 ? 4 

on  t i ?o 

0 0?  s 

10 

nr  2 1=1,6 

IT"  002  6 

Y =°  ( I ) 

£ 0027 

F < I) =Z  3(  Y.C'fcl  1 

t 002  0 

2 

F( I )=Z4(Y,DEL ) 

“ 0029 

a i =r  = i.. 

007  0 

A 2=  1 . - 3 “ L 

00^1 

20 

T R L = 0 . 

jT~ 0032 

90  3 I =1  ,3 

l 003  3 

xi=c(2*i-n 

l 0034 
003  5 

X2=E(?*I ) 

Y 1 =P ( 2*1-1) 

0036 
003  T 

Y2=f (2*  n 

C = A 1 = ( RT-.  ( AR,  XO  ,Y  J, 
C = A 7# « F7°  ( AP, XO  0, 

JT  0038 

i 0039 

3 

TR  L = T P L +W ( I 1*  (C+0 ) 

0040 

RE  TURN 

0041 

ENr 

10 


000  3 
000^ 
000  5 
0006 
0007 
0008' 


c J - ' < ‘ ’ ’ I ' V L I [ 7 T f r 7 4 1 4 < ; 2 1/55 


00  J 1 
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0010 
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■i  > 1 1 
001)2 
0003 
000  A 
0 00  5 
JOCK 

0 JU  7 
OOOH 
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001  1 
0 '1  2 

001  3 
00  1 A 
0015 

00  1 5 
0017 

001  8 
0019 
00?  0 

002  1 

no?? 

00  2 3 
u J 2 A 

00Z5 

0026 

00?  7 
0 0?  8 
00?  9 

17030 
003  l 
0032 
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0035 
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1 1-1  , T X 1 1 t /i ) 


Cl  12) , TAU( 12 , 2)  ,PAT ( ?0 ,1 3 ) , ANS ( 1?  ) , A ( 


AND  SPANW1SE 


,0.  ) 


S' HP  1J1  I ■;  T •,nc  ( ^ i ; , 

DIMENSION  11(20) 

DIMENSION  C01K9  ,A0)  , 

REAL  MU 1 , MU  ? 

PP  ?0  1=1,12 
A { 1 ) =0  . 

TALK  I , 1 ) =0. 

?0  T A U ( I ,21=0. 

FIRST  Mpnc  CHORDw! SE 
N=NSD-1 

P AT(  ] , l >=  1 .-rpij(  1 , 1 I /rr(  1 ) 

P AT ( 1 ,2)  =1  . 

P Mil  , 3)  =1  . -c  CU(  1 ,2  l/CC  ( 1 ) 

no  1 I = 2 , N 

RAT! I , 1 ) = 1 - -COU  ( 1 , I )/CC( 1 ) 

F AT( I , 2) =1.-C0U« 1, I- 1 ) /C0< 1) 

1 RAT(I,3)=1.-C  DU (1, 1*1 ) / CO ( 1 ) 

CALL  r-LSO(R  AT  ,ANS,IL,N,2,BUG,0 
Al  = . 5* ANSI  l ) 

SQ=SQRT ( Al^Al+ANS (2) ) 

M*J1  =A  1 •‘•SQ 
MU2  = A 1-5Q 

TAU(1,1)=ALDG(MU1  l/H 
T A U ( 1, 2) =ALOG(MU?»/H 
DO  2 I = 1 , fvc  0 

2 A 11 ) =A  ( 1)  M ) .-CCU(  1 , I )/CH  (1  ) -MU  1**1  ) / ( MU 2**  I -ML  1 **1  ) 
A(  1) =A( 1 J/NSD 

SFCOND  MODE  CHORDW  I SE 
005  II =1 ,NWP 
J=NMS  + I ? 

TAT-l  . 1 ) r 1 . - . 5*  COl  I ' J , 2 ) /C C ( 

RAT\  1,21=1. 

RAT (1,  31=  l.-C0U(J»2)/C0(J) 

DO  3 I =2  , . 

path  , 1 1 = ,-rnur  j , M /cc(  ji 
p AT(  ! , ? )-- 1 .-mu  ( j,  t - 1 > /cn;  j ) 

3 p at ( i ,3) =1 . -c ru( j , i +i ) /cr ( ji 
RAT(2,?)=PAT<  1,1) 

CALL  r | s 0 < P at  , A NS , I L , N , 2 , BUG 
A 1 5*ANS(  1 ) 

S Q = S OF  T ( A 1 *{■  1 •‘•AMS  (?)  1 
Mill  = A It  CO 
Ml|?=Al-c° 

TAUIJ,  1 1 = ALPG(MU1  ) / H 
TAU( J , 2 1 =AL  CG ( MU2 ) /H 

A(  .1)  =(  1 .5*fOU  ( J , 2)  /C0(  J)-MU1>  /(  MU  2-MU  1 ) 

EC  A I =?  , N c n 

A A(.J)  -{.  ( .!)♦(  l.-Crn(J,  I )/CO(J)-vl  lft’  1 ) / ( MU?*f  I — A L 1 
5 MJIU.IJI/'";^ 

RETURN 
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SAIL-RUDDER  INTERACTION  RESPONSE  PROGRAM 
INPUT:  1)  Vehicle  Geometry 

2.)  Trailing  Vortex  Sheet  Angle 
3)  Exponential  Approximations  to  Circulation 
Response  of  Sail 

OUTPUT:  1)  Step  by  Step  Forces  and  Mo  ments  on  Rudder 

due  to  Unsteady  Sail  Wake 
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. * ■>. 
•V> 


V r 1 


2 1 


MA  I N 


MAT  r 


1 h 2 6 ( > 


1 U < 


C I vrNS  I r N XI  A j)  .PTUill  , CP( 20, 2)  ,P  VP (20,2)  ,0W  I ( 40,  40  , 3 ! 


n I '•KNS  If  N I’wR  (40,40,  3 ) , DWR  1 < 20,  13)  ,0WR  ( 20, 1 3 ) ,R(  20)  ,C0l.  ! 1 2 , 4 . t 
DIMFNSICN  SF ( 5)  , R M { b ) 

I)  I Mf  NS  II  N P'JF  ( 40)  ,PIJM(40)  ,RUR  (40)  ,FM(2,40),CFM(3,40),DFM(3,40( 

0  I MFNS I GN  I L I 2 0 ) 

RPA01  00  ,NMC  ,NMS  ,NCPC  , NCPS,  NSD,NX,  TH,RHR  ,PRD,RL  ,HCP  ,RHS,  RSL  , PCJW  , , U 
1RN 

100  FORMAT ( 61 5 » 9F  5 . 3 ) 

REA01 LO,FVS , A , T1 , T2 
110  FORMAT (4F12.4) 

NPTS=NCPC*NCDS 


N MO  S = N M S = N MC 


REA0200, CCP<  I ,1)  ,1  = 1 ,NPTS) 

200  FORMAT ( 16F5. 3 ) 

ARR=.5* (RRD-RHP)*(RRD*RHR) /RRD/HCR 
Yl=l.-(RRC-ftHR)/ARR/HCR 
ARS=.5*  (RSL-PHS)MRSU-RHS)/RSL 
YH=.6366  1 9 7 * A R S I N ( Y 1 ) 

00  1 I = 1 , NC  PC 

1 CP( 1M I-l)*NCPS,2)=(Yl+. l)*ARR 
DO  2 1 = 1 , NC  PC 

no  2 J = 2 , NC  P S 

2 CPIJM  1-1 ) r NC  P S « 2 l=(2*J-3)*.5*ARR/(  MCPS-1 ) 
ROR=RRC-APP*HCR 

R OS  = R SI -ARS 

00  3 I = 1 , NP T S 

R VP ( I, 1)=CP( I ,1)*HCR+RL 

3 RVP( I ,2 ) =CP( 1 ,2) -HCR+ROR-ROS 


X ('.)  = 1 . 


9 T ( 1 ) — P S L 


CAIL  WCTIBCV,, STERN, X, RT  ,NSD> 

025  FORMAT ( ‘O' , 2X  , E 12. 4, 2X , E 12 . 4) 

CO  4 I = 1 , NP  T S 
XO=RVP<  I , 1 ( 

Y 0=  R VP ( 1,2) 

CALL  UTRLD( XO ,YO,TH, X,ARS,RT , l ,NX ,NSD, I ,DWI ,RHS, RSL ) 

4 CALI  UPN0< XO, YO, TH-X.RT, ARS, 1 ,NX,NSD, l ,DWB,RHS,RSL ) 
DO  5 1=1 ,NPTS 

PR  I NT  300,  (DWI  ( I , J, 1 ) , J= 1 ,NX ) 

5 PR  1 NT  300  , ( D W 13  ( I , J , l ) , J=  1 , NX  J 
CO  6 I=1,NPTS 

XO=CP( 1,1) 

YO=CP( 1,2) 

CO  o I J = 1 , N M C 


00  6 J = 1 ,NMS 

DWPI(:,(!J-1)”NMS+J)=SINT(X0,Y0,J,IJ,ARR,YH) 

IF!  IJ.EQ.1 ) CWP 1 ( I , J)=0WR1  ( I , J )+TRL ( ARR, YH , J,X0,Y0, 0, ) 
( C ( N T I N U F 
PRINT  -0 

9 )0  FORMAT ( ‘ C'  , 2X  STATEMENT  1*1 


cr  o 
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Pc  I X ( 1 )»PT(I),I=1,N$D) 

Or  7 I = i ,NPTS 

DWR 1 ( I , NMD  S + l ) =0^8 (1*1,1) 

NE=NM0S+1 
rn  7 j = i , nf 

CWR(  I,  J)=()WR1  ( I , J ) 

C ALL  01  SC ( DWR  * B * ! L , NPT S , NMDS , PUG , 0 „ , J.  ) 

DO  8 1=1, NMDS 
C0U( 1,11=0(1) 

00  6 1=2, NX 
OCJ  10  J = 1,NPTS 

DWR  1 ( J , NMUS  + l ) =DWR  1 ( J,NwDS+iM-DW0(J,i,l)-CWB(J,  1-1,1)  + CW  I ( J , 1 , 1 ) 
OWP ( J , NMCS+ l ) = OWP l ( J , NNDS  + 1 ) 

DC  10  I J = 1 , NMDS 
DWR(J,IJI=DWR1(J»IJ) 

PRINT  400, I 

FORMAT ( 'C' , 5X , •DOWNWASHFS  AT  STEP*, 13) 

PR  1 NT  300, (CWR ( IK ,N^DS+1 ) , I K = 1 ,NPTS) 

CALL  GLSQ(0WR,B»1 L , NPT  S , NMDS , 0UG , 0 . ,0. ) 

00  9 IJ=1,NMDS 
COU(  IJ,  I )=B( I J) 

PR  INT950 

FORMAT! 'O' ,2X, 'STATEMENT  2') 

PR INT925, ( X( I ) ,RT( I ) , I = 1 ,NX) 

PR  INT350 

FORMAT ( • O' , 2X , *C I RCULAT I CN  RESPONSE  OF  RUDDER  TO  UNIT  STEP') 

CO  11  1=1, NMDS 
PR  INT375, I 

FORMAT ( « 0 • , 2 X , ’ MODE « , I 2 ) 

PRINT  300 , (CCU( I , J)  , J = 1 , NX ) 

FORMAT! • C' ,2X , 10E 12. 4) 

CALCULATE  FORCES  r.  MOMENTS 
CALL  SR( ARR,YH,NMS,SF,RM) 

DC  12  ! = l » N X 

RUF(  I )=-COU(  1 , I )*3. 1416*2. “SF  (1  ) * HCR “ + 2 
RUW<  I ) = 2 « * S F ( 1 ) * 1 , 57  08*COU(  1 , I ) * HC  R * * 3 
PUC(  I ) =-COU ( 1 , I)*3. 1416*2. *RM(1 ) * HC  R*  * 3 
CO  13  J=2,NMS 

RUF  ( I ) =RUF  ( I ) -C.OU  <J,I)*3.1416*2.*SF(J)*HCR**2 
RUM ( I ) =RUM( I ) +2.*SF ( J l*l.5708*C0U( J , I >*HCR**3 
R UR  ( I ) = RU°  ( I )-3.14  16-t2.*RM(  J)*COU(  J , I ) * HCR*  *3 
SECCNC  MODE  CHOROrfISE 
CO  12  J = l , NMS 

RUM(  I ) = RUM(  I )-.  7864»C0U(  J+NMS,  I (»2.*SF(  J)  •»HCR**3 

DO  18  1=1, NX 

R UM ( I ) =RLM(  I l-RUF (I)  »RL 

‘-UP  ( 1 ) =RUR  ( I ) +PIJF  ( I ) M RRD-ARR*HCP  ) 

PR  I NT  380 

FORMAT (• O'  , 2X  , ‘ FORCF  RESPONSE  CF  RUDDER  TO  UNIT  STEP') 

PRINT  39  0, ( X ( J ) ,FUF ( J ) ,RUM( J) ,RUR( J ) , J=l ,NX) 


^ tnnt'.Tt..-  ■f-  w.fj- 


, . «»!*?**?«*■  i . 


C \< 
< . 'i 


'>«’  IV  v'  I : V F l 2 1 


MAIN 


0 A T £ = 74260  17/i 


9 
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2 

A 

5 

6 
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S 
9 
0 
1 
2 

1 

A 


>9  0 F i PMAT!'0',2X,'X  = *,E12.4,2X,,SF  = ',E12.4,2X,I’VM=«,H12.4,2X,,RM  ’ 

12.  A ( 

CO  15  1 = 1,  NX 
FM(  1,  I ) = RUF  ( I ) 

FM(  2,  1 ) =RlJM(  I I 

15  FM( 3,  1 ) =RU9  l I ) 

H = X ( 6 ) - X ( 5) 

no  16  1=1,3 

CALL  Of RY(NX,HfFMfDFM,  1) 

CALL  CONV(NX,H,DFH,FMfA, T1,T2,CFM,  I ) 

DO  16  J = 1 , N X 

16  CFM( I ,J)=CFW( ! , J I * 3. lA16ttFVS 
PR INT395 

395  FORMAT C O' , 2X , 'CONVOLVED  FORCE  L MOMENT  RESPONSE'* 

PRINTA10,  ( X ( I *,CFM(l,l*  ,CFM(2,1  ) , C F M ( 3 » I * ,1=1, NX* 

A10  FORMAT ( 'O'  ,2X , ' X = * ,F 1 2 . 4 » 2X SF  = ' ,E12.A ,2X, ‘ YM=  ' ,2X,F12.A,2X»  *RM; 
1.F12.4) 

STOP 

END 
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FuNCT ITS  S i NT  ( X , Y t N » M , A R , V H I 
r I men SIGN  R ( 6 ) , E I 6 ) , w( 3 ) ,F(6),CL(6)  , C 2 ( 6 ) 
Z 1 ( A, 0 ) = 3‘  ( l . + A ) -1 . 

1 2 I A » 3 ) ’-  B r ( l • - A ) ♦ A 
Z3(Af6)=A*(l.-M3)+0 
Z 4 < A,3)=R*  ( 1 . - A ) +•  1 . 

ZETA  = . 6 366197 'AOSINIY/AR) 

R (lt=. 2386192 
o(2I=-l  .AR(1> 
o( 3)=. 6612099 
c 16 )=-l . *R I 3 ) 

R I 5 ) = . 9326695 

R I 6 » =- 1 . *R (51 

W I 1)  = . 66  791  39 

W ( 2 ) = . 3607616 

W (?)=.  1 713265 

IF( ZETA.LT.O. 1G0  TO  3 

C1  = ZETA 

02=1. -ZcTA 

DO  1 1= 1,6 

6 1 = R ( I ) 

E( I l = Zlt A1 ,ZETA) 

E ( I ) = Z2( A1  , ZE 1 A) 

8 1 - F ( I ) 
p 2=  r 1 1 1 

E TA  1 = AR  *SIN(  1 . 5708*81  ) -Y 
ETA2-AR*SIN( i« 5708*82 )-Y 
C1I  I)=C  INTIX.ETAl.MI 
C2 ( I I = C INTI  X , E TA2  ,M> 

GO  TO  5 
ci=i.+zsta 
r?--l -- ZF Ta 

DO  6 1—1,6 
A 1 = R I ! ) 

I <I)=Z3(AL,ZFTA) 
f I I l = Z6IAl,Zf  T A I 
3 1 - E I I) 

F)  2 = F I I ) 

ETA  1 = AR  * S IM  1 . 5708*31  )-Y 
FTA2=Afi*SIN( 1.5708*B2(-Y 
C 1 (I)  =C  INT ( X, E TA1 ,M) 
r? i n = r inti x , e t a? ,m) 

G 1 = C. 

G2  = 0. 

no  2 11=1,3 
J J = 2*  I I 
( J =2*  I 1-1 
X 1 = E ( IJ  ) 

X 2 = E ( J J ) 

Y 1 = F ( [ J ) 


■ ■ 


— ,,i'iH^®'J?fRtt!3*«tTnw/(wtv'i«.i«a"»i 


t.:G< 


D AT  F = i‘i? 6 6 


Y2-F  ( 
1 = C 1 
W2  = Cl 
R 1 = C2 
R?  = f.2 
G 1 = C 1 
G2  = G2 
S ! N T= 
RE  TUR 
END 


\ \j  i 

< I J ) 

( J J > 

♦ Will  )M  FI  ( XL  ,Y,  ml  , AH  ,N,  YH»*FK  X2,  Y,  W2  ,AP  ,N,  YH)J 
«-W{II)'»(El(Yl,Y,Rl,AR,N,YH)+Fl(Y2,YfR2,AP.,N,YHn 
(G1*01+G2*D2 )*. 19635 

M 
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10 


20 


HJMCT I C \ C!NT(X,FTA,M) 

rP'ENSKN  R(6),t(6),F(6),M3) 

Z 1 ( A,B ) =A- ( 1. *R)  + R 
Z2(  A i 8 ) = 0*  I l . - A ) <■  1 . 
Z3(A,B)=R-(1.+A)-1. 

Z A ( A , R )=BM  l.-Al+A 
L F l =, 6166197 “ABSIN(X) 

R ( 1>=. 2386192 
P ( 2 1 =-  l • “R ( 1 ) 

R<3)=. 6612094 
R (4J=-1  . *R( 3) 

R ( b 1 = .0324695 
P (6)  =-l  . *P.(  5) 

W(  1)=:.  46791  39 

M2  ) =.  36C7616 

W ( 3 ) = • 1 7 1 3245 

I F ( DFL • GT . 0. ) GO  TO  10 

HO  1 1=1,6 

Y=R ( 1 ) 

F ( I ) =Z 1 ( Y »DEl  ) 

F ( I ) = Z 2 ( Y , DEL  I 
A 1= 1. +DEL 
A 2 = - 1 . * DEL 
GO  TO  20 
DO  2 1=1,6 
Y=R (I) 

E ( I )=Z3< Y.DFL ) 

F 1 1 ) =Z4 1 Y ,OEL  ) 

A 1 =DE  L 
A 2=  l.-OEL 
C I N T = 0 , 

DO  3 1=1,3 
X L = E ( 2 * [-1  ) 

X 2 = E ( 2 r 1 ) 

Y l = F ( 2 * 1-1 ) 

Y2  = F ( 2W  l ) 

D = A1M  FC(X1 ,ETA,X,M)+FC(X2*ETA,X,H) ) 
C=A2*(FC(Y1,FTA,X,H)4FC(Y2,ETA,X,M11 
CINT=C  INT+M  I 1 * ( D+C  ) 

R F TUB  N 
END 


fcS- 


F W [ V 
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fc 


U A T F = 7 9 2 bo 
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> I 


4 

5 

7 

8 
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TUNC! I ON  FC  (X  ,F1  A,XO,MJ 
IF (M.EO.O)G=SIN( l.5708*X) 

I FI M.FU.I IG  = 1 .-SIN(l .5708*X) 

I F { w.FO.2 )G=- .5* ICGSI 3, 1 4 1 6” X )+SIN( 1 . 57  08*X ) ) 
IFIM.GT  ,2»r,=  .5*<  COSCM*l  .5708*  I l.+XM  ^CCSI  (M-?  1*1. 
A = S I N I I. 5 70 8* X) 

FC  = G*SORT ( ( A-XO I * ( A-XOI  +ETA*ETA  )/ ( A-XO) 

RETURN 

END 


5708*  I 1 + X ) ) ) 


?s 


a ^ ^ 


OAT  F 
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FUNCT  I ( N F l ( A,Y,8,AP,N,  YH' 

IFIN.FG.  1 » Gil  TO  1 
IF(N.E0.2)f»C  TO  2 
I F ( N.  E 0 . 3 I GO  TO  3 
CONTINUE 

IF(A.LT.YH|C-S=C0SU.6708“(AM.>/(YK+1.  I ) / (YH+1  . » 

IF(  A.Gl  .YH)GS  = CPS(  1. 5708-  1 A-l  . ) / ( YH-l.  »)  / < YH--1.  ) 

GO  TO  4 
CPN1 INUF 

IF(A.LT.YH)GS=Ci!S(’ .5708*(A+-1.)/(YH+1.))/(YH+1.) 

I F ( A.CC.Y8  >GS  = COS( / . 5708*A/YH)/YH 
IF1A.G-.  .0.  IGS  = -2.  * CCS  ( A- 3. 1416) 

GO  TP  4 
CONTINUE 

tF(A.’T.yH)GS  = 3.'*CnS(4.7124*<A+l.>  / I YH+  1.  ) ) /<  YH+1.  ) 

IF( A.Gi  . Yu } GS  = 3 • - COS (4.7124s ( A- 1 . ) / ( YH-l. ) ) / I YH- 1 . ) 

F1  = H»-GS/(AP*S1N(AM.5708)-*I 

return 

ENO 


5fO< 


P AN  [V  c,  l ! V F I ^ 1 


LATE  = 76  26^ 


7/  Jy IP 


F lli'JC  T I f N TPL ( AR , YH, J , XO, YO , TH) 

DIMENSION  R < 6 I ,E  (6)  , K6  I ,H(  31 

Zl(4,3hA‘ll.«R|tB 

Z21 A,  8 I =9" ( I. -A ) + l. 

7 3 ( A,  B ) = 0’'  ( 1 . + A ) — 1 . 

Z4(  A,  8)=0M  l.-AI+A 
rEL=.63fcfcl97»AF'SIN(Y0/AR) 

R<1>=. 2386102 
WI2)=-1.»R( 1» 
p ( 3)  = .661209A 
R 1 4)=-l .*R( 31 
R<51=, 9324695 
R(6)=-1.»R(5) 

W(1  I = • 46  79 1 39 
W ( 2 ' = « 36C 7616 
Vv  I 3 I = . 1713245 
TF1CEL.GT.C.1GO  TQ  10 
CO  1 1=1,6 

Y = R 1 I ) 

E 1 I) =Z 1 ( Y » DEL  I 
l F ( I » =Z2  ( Y »DEI.  5 
A 1=  1 . «-DEL 
A2=-l.*0El 
GC  TO  20 
10  CQ  2 1=1,6 
Y=R ( I } 

El  I l = Z3 ( Y ,OEL ) 

2.  F<  1 1=Z4(  Y.DEL  ) 

A 1=0EL 
A2=  l • -DEL 
20  TRL=0 j 

DO  3 1=1,3 
X l = E( 2* 1-11 
X2=E ( 2* I » 

Y 1 = F ( 2r  I-  1 1 
Y2  = R( 2®  I ) 

R = A1* I FTRI AR, XO,YO,YH, J,TH,X1  H-FTR1 AR,XO,YO, YH, J,TH,X2) I 
C=A2* ( ETRl AR , XO, YO.YH, J, TH,Yl 1 +FTR < AR , XO, YO, YH, J , TH, Y2 ) ) 
3 TRL  = TPH-l,< I )*  <C  + D) 

RETURN 
1 NO 
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HJNCT  luN  tSPAM  A.N.YHl 
I r (N'.fcO  . i )CC  TO  1 
I F{  N.  t O.  2 I 00  TO  2 
!F|H.EQ.?!(in  TP  3 

1 continue 

P ( A.LT  .YH)G.S=C0S(1  . 5 708*  ( a + 1 . ) / ( YH*1.  ) )/  ( YH*1.  ) 

I F ( A.C.F  .YH)GS=CPS(  I.S708M  4-1  . ) /(  YH— 1.  )>/  ( YH-1 . > 

GO  TO  4 

2 CONTINUE 

If  1 A.LT.  YH)GS--cOS(  U 570a*<  A+l.  ) / { YH+I.  ) )/ ( YH+1.  I 
I F( A.GE .YH)GS=COS( 1. 8 7 03* A / Y H ) / YH 
I.'l  A.GF  .C.  ) GS  =-2  . *CDS  ( A*  3.  1416) 

GO  TO  4 

3 CONTINUE 

IF( A.LT. YH)GS=3.”C0S(4. 7124*1 A+l. )/ ( YH+1. ) )/l YK+1  . ) 
IFlA.Gt. YH)GS=3.#C0S<4. 7 1 24* ( A- L , ) / ( Y H- 1 . ) ) / ( YH- L . ) 

4 FSPAN  = GS 
RETURN 
END 
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35 


32 


30 


51 


50 


52 


GLSQ 
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SURRTJTlNfc  GLSQ(A,Xf IL,N,M .ALPHA, cl»E2) 
DIMENSION  A ( 20  * 1 3 ) • X < 20  ) . IL(20i 
M H = M •*■  1 
L L - 1 

00  60  J=l.MM 
U ( J)=0 

1 = 1 

CO  3 K=1,MM 
1 1 = 1 + 1 
DO  A J=  1 I * N 

IKABSUCJ.K)  )-El)4,4,6 
t 1=  SORT  ( mJ,Kll:*;,2+(AII,KII**2) 

S = A I J , K ) / T l 
C=A! I ,K ) / T 1 
CO  5 L=K,MM 
T2  = C*A(  I ,D+S*A(  J,L) 

A ( J , L ) = -S’>Al  I ,l)+C*AIJ,L) 

A ( 1 * L I = T 2 
LL=LL ♦ 1 
CONTINUE 

I F ( ABS ( A( I ,K1  » — E 2 >3, 3,8 
I L ( K I = I 
1 = 1 + 1 
continue 

X(MM) =-1.0 
II=M 

00  35  1=1. H 
X(  11  = 0.0 

DC)  30  J = 1 1 M 

1 F ( I L ( 1 1 1)30, 30, 31 
S = 0.0 

LL=  1 1 + 1 
I = ! L ( I I 1 
CO  32  K = LL 
S = S + A I I,K5  + X(K1 
xmi=-s/A(  i, in 
H = I l-l 

IF(  IL I MM ) 150,51,50 

ALPHA=0.0 

GO  TO  52 

I = I L ( 1 

ALPHA  = A ( I . MM  1 

RETURN 

END 
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SUP -HUl  T 1 A fc  U^NO'  XD,  YO,  T , X,RT  , AP  ,NMS>  ,NX  ,NSD,  I C P , DWB  , RM  S , R SL  ) 

0 I^FNS ION  X (40  I . RT (431  , DWB( 40 ,40,  3) ,R( 6) , E( 6) ,F(6 ) , W( 3 ) 

7 1 ( A,  B,C  ) = .5*  ( AM  C-B  ) +C  + 8) 

H = X ( 2 ) -X (1) 

Z0=(X0-1.  » - S I N ( T > 

xoi  = cxo-i.  >*r.osm 

CO  1 l = 1 , NM  S 
CO  1 N = L , NX 
I F( N.GT .NSD1GO  TO  6 
XA=  X < N ) - 1 . 

PTE=RT ( N) 

XH=  X ( N ) 

P H = HR  A D ( X H ) 

ARN=AR*(RTE-PH)/(RSL-RHS) 

Y0L=Y0-RSL+AR4-RT£-ARN 
YHl=AP.  N-RTE  + RH 
YHT=ARS INI YH1 /ARN )♦. 6366197 
6 CONTINUE 

IF(N.GT.NSD)X(N)=X(N-1»  *•  H 
I F(  N.GT  .NSD  I XN  = XN«-H 
I F ( YOl.GE.YHl ) GO  TO  2 
HI  = .1*(  l.-YHT) 

Y 1=  YHT  + H I 

VINT=FU3ND(YHT  ,X0l,Y01,Z0,XN,APN,YHT  ,1  ) <-4  . •*  FIJBNDI  Y1  ,X01 , YC1  ,Z0, XN, 
LARN.YHT, I ) + FU8N0{ 1.  ,X01 , YOl.ZO, XN, ARN ,YHT , I ) 

00  3 J=2,8,2 

Y 1 = YHT  + J*HI 
Y2=YUHI 

3 VINT  = VJNTf2.‘FUBND(Yl,X01,Y01 , Z 0 , XN, ARN , YHT , I ) +4 . *FUBND ( Y2 , XO 1 , YO 1 
1,Z0,XN,mPN,YHT,  I J 

DWB ( I Cp , N , I ) =Hl*VINl/3. 

GO  TO  L 

2 CIF=.6366197*ARSIN(Y01 J-YHT 
UP=  YHT+  2 . *D I F 

1FIUP.GT.1. IUP=2.+(DIF+YHT)-1. 

C>  L=  ( UP-YHT  I * . 5 
02= ( l . -UP ) * . 5 
R(l)=. 2386192 
R(2)=-l.*R(l> 

R ( 3 ) = . 6612094 
PC4)=-1.*R(3) 

P< 5)  =.9 3 2469 5 
R *6>  =-  I..«R(  5) 

W(l)=. 46791 39 
W<2>=.  3 6 C 7 6 16 
W ( 7s  I = . I S 13  245 
00  4 i I = 1 ,6 
A 1 = R i IT) 

f 1 l 1 5 = 2 1 ( A 1 , YHT, UP) 

4 F m=Zl(  A1  ,UP,1.  ) 
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FUNC  Tit  N FUBNO  ( V , X'J,  YO,  ZO,  XN,  AH  , YH,  I ! ) 

IF  ( I I . FQ.2  >GO  Ttl  1 

I F ( II . FQ. 1 ) GS=S IN ( 1 . 6 708*  < Y-l . ) / ( YH- 1. ) ) 

I F ( II  .FQ.3)0S  = SIN(4.  7124MY-1.  )/(YH-l.  ) ) 

cr  TO  2 

IF( Y.GT.YH) GS=S INI  1.5708* Y/YH) 

!«;(Y.GF.C.)GS  = -S1N(Y*3.1<»16) 

CONTINUE 

ET4=AR*SIN(Y51.5708)-Y0 

CEN  = FTA>ETA  + /0<,Z0+{XN-X0)MXN-X0) 

FIJRN!)=.  L 2 AF  • ( XN-XO  ) » G S*CPS  ( L.  5708*  Y ) /OE  N*  * 1 
RE  TURN 
f NO 


IV 


— 


<>7< 
•_  > 6 


l b V-  I <r  1 


I) T ;? l T 


HATE  = 7426o 


17/33/18 


hUtMGt'T  I\r  IIT  RL0  1 XO,  YO,  1 , X,  A8  ,«  T .N^S  ,NX  ,NSO.  1CP»DWI  ,RHS  ,P  SL  ) 
r-  I m F m s It  N X {6  J)  , *3  T < <t  0 1 , CM  ( 40  ,40,  5 ) , R ( 6 I , F ( 6 ) , F ( 6 ) , W ( 3 ) 

Z 1 ( At  <},(.  ) = .5M  A«  < C-8)  +C  + « ) 

PRI\TLCJ,XC,YO,T,AR,NN'SyNX,NSD,ICP,RhS,RSL 
Ft  R MA  f < ' 0 ' » 4C  12.4,4l5,2F12.4) 

H=X ( 2 ) -X ( 1 ) 

2 C= ( X0- 1 . »<  SIM 7 ) 
x o i = ( xo-i . ' * f. o s ( T ; 

DO  1 I = l , N M S 
no  1 N = 2 < NX 
I F ( N . G T , N S 0 ) G i J TC  6 
X8=  X ( N ) - 1 . 

XC  = X( N1-.5*  H 
X T = XN-H 

P TC  = ( ° T (M  + RT  (N-l  ) )»  .5 
PHC-HR AC( XC ) 

ARN  = A-tM(RTC-RHC)/(RSL-RHS)l 
Y01  =YO-RSl+AR+RT  C - A R K 

yhi=ar\-rtc+«hc 

YFIT  = AR  S IN  I YHl/ AON)*.  0366  197 

> contimue 

IF ( N.GT  .8 SO )X  (N ) = X ( N- 1 ) +H 
I F ( N.GT  . N SO  ) XN=XN*H 
I F( N.G T .NSD ) X T = X T f H 
IF  I YOl.GE.Yhl >G0  TO  2 
H I = . 1 - ( l.-YFT  ) 

V1=YHT+HI 

Y INT  = FUTRLD  (YHT  ,X01  , Y01  ,ZO,XN?XT,  ARN,YHT,  [ ) * 4 . * FUTR  LD  ( Y 1 ,XOl,Y01 ,Z 

10. XN,  XT  ,ARN,YHT  , ! M-FOTRLDI  1 . , XOl.YOU  ZO  ,XN,XT,  ARN,  YHT,  I ) 
rn  3 J = 2 y 8 , 2 

Y 1 = YHT  + J"H I 

Y 2 = Y 1 +H I 

VIMr  = VINTf2.»FUTRLD(Yl,X01,Y01,  ZO  , XN , XT  , ARN  , YFIT  , I)  +4  . * FUTRLO  ( Y2  , XO 

11, Y01,Z0,XN,XT  yAPN*YHT , I ) 

CWI  I ICP.N,  I )=HI*VlNT/3. 

GO  TO  1 

! 0If-=.6366l97*'ARSIN(Y0l)-YHT 
UP=YHT+2.*D IF 

IHUP.GT.l.  )UP=2.*(DIF  + YHT»-1  . 

01=  (UP-YFT  ) * . 5 
02= ( 1 . -UP  I “ .5 
R ( 1 ) = . 2386192 
R(2)=-1.*R< 1) 

R < 3)  = . 6612094 
R(4)=-1.*R(3) 

R ( 5 J = .9324695 
R ( 6 ) =- 1 • aR ( 5 ) 

W ( l >=,4679  L 39 
Vy  ( 2 I = . 36076  16 
W(?  ) = . 1712  245 


‘At  Hit'  <Ztl  i 


23< 


I 


15  AM  IV  C.  L f V r L 2 1 


iJTPL  D 


DATE  = 7A266 


1 7/ 3 ^ / l 8 


C*1  A 11=1,6 
A 1 = R ( I I ) 

t (II)  = Zi(Al,YHT,IJP> 
r I I I ) =Z  n A 1 ,UP,  1 . ) 

r i = o . 

G 2=3. 

rn  5 II  = 1 , a 
X1=E(2*  I 1-1  ) 

X2  = p I 2*  I I ) 

V.  1=M  2-  ! 1-1  ) 

W2=F<  ?-■  I I ) 

G1  = G1*W(  II)M  FUTPLIM  XI  ,X01  , Y01  ,ZO,XN,  XT  , ARM,  YHT  , I)  i-'-'J  T PI  0 < X 2 , X0 1 ,Y 
10 1 , ZO , XN, XT ,ARN,YHT, I ) ) 

G2  = G2*WIII IM  FUTPLDI W1 ,X01 , YO  l.ZO.XN, XT ,ARN ,YHT, I ) +FUTPLD  (W2iX01 ,Y 
101,Z0,XN»XT ,ARN,YHT, I I ) 
own  I CP  ,N,  I)=  (Gl*Dl  + G2*02) 

CCNT INUE 

RETURN 

END 


R AN  IV  1 1 t.  r V rL  2 1 


"UTPLO 


DATE 


74266 


L7/33/1H 


1 

2 

3 

4 
6 

0 
7 
H 
9 
0 

1 
2 

3 

4 


FUNCTION  FUTRt_9(  Y,XO,  YO,  ZO,  XN,XT  , AR  , YHf  II) 

IF( II . EQ.2IGU  TH  1 

IM  II  . CQ . I )GS=COS ( 1 .6708* I Y-L . ) /{ YH-1  . ) )/ l YH-i  . ) 

IF(  I I . FQ.  3)  GS  = 3.’lC(JS(4 . 7 124M  Y- 1 . ) / I YH-1 . ) ) / ( YH-l  . ) 
GO  TJ  2 

1 IF < Y. GT.YH )GS=COS(  1.5708  ‘Y/YH1/YH 
I F ( Y.GF.  .0.  ) GS=-2.  “CO  SI  Y’  3.  141  o) 

2 CONTINUE 

ETA=AR“SIN(Y'-1.57  08)-Y0 
DEN=ETA3FTA+Z0*Z0 

F l=( XN-XO )/SCRT< I XN-XO) » (XN-XO) +DENI 
F?=(XT-XO)/S(JftT(  ( X T - X 0 ) * ( XT-XO) +DEN  ) 

FUTRLD  = -. I25-GE»  ETA* (F1-F2) /OFN 
RETURN 


^ . 
• ' - A. 


DATE  = 7426c 


1 7/ 3 3/  1 8 


SLP-^U'ITlKr  wCT ( A , 3 , X , R , N ! 
fJI^FNS  I IN  X(*0)  » R ( 4 0 ) 

F'NT(EtFfG,I,0)  = E/  ((r-l.-<I-l>*niMF-l.-<I-l)*D>*G’*G)**l.5 
h= .025' ( E-A ) 

HT= (R- 1 . 1 /N 
N ! =N+  1 
CO  1 J = 1 ,N1 
R0=  R ( J I 
Fl=ASLP(A) 

A 1 = A-*-H 

E2  = ASI.P  ( At  ) 

F 3=  ASLP ( R I 

X!NT  = FINr{Fl,A,R0,J,HT)^4.-FINT(F2,Al,R0,J,riT)rrINT(F3,B,R0,JtHT) 
CO  2 I J r2  * 38  * 2 
A 1 = A+ I J * H 
A 2=  A 1 + H 
F1=ASLP( Al ) 

F2=ASLP ( A2 ) 

XINT=XI\T+2.“FINT(E1,A1,R0,J,FT)+4.*FINT(E2,A2,R0,J,HT) 

X INT=H'HT»X INT/18.85 
P<  J-*-l)=R(J)M  l.+XINTI 
X{ J+l) - X ( J ) »HT 
RETURN 
FND 


p am  iv  r,  i.  fvel  ?i 


ASLP 


DATE 


7 A 2 6 6 


1 T/33/18 


FUNCTION  ASLP(X) 

BPW=-5. 125 
STERN=7.031 
X l=-4. 25 
X 2 = . 3 7 5 
X3  = £. 

R 1 = .4  3 5 1 
R2=.8438 
R3=. 328 l 

IFIX.LE.XIIGO  TO  1 

1 F ( x .GC  . X 3 ) GO  Tn  2 

A=<  Ri*R  3-2.  *R2)  ".5/ (X2-X1)*  *2 

B=( R1-P3I/IX1-X3) 

C = R2 
XD=X-X2 

ASLP  = 6.  2832*(2.*A*A*XD»XD«XD+3.  *■  A*  P*  XD*  XD+ 1 R*B+2r  *A*C  ) * XD+C*B> 
GO  TO  3 

1 ASLP=3. IA16*R1*R1/(XL-80H) 

GO  TO  3 

2 A SL  P = 3 . L416*R3+R3/( X3-STERN) 

3 CONTINUE 
RETURN 
END 


4 

* > * 


l f V-  l 


= 74266 


17/33/18 


SlI-AfUir  1 Nr  OF  ' Y ( N'  ,Hf  f-  , OF  , 1 I ) 

DPIFNS  If  N f ( 3 .40)  ,01-  ( 3,40) 

Fl  = 1 - / H 

nf  ( ! I ,M=H|  * ( P ( ! I ,N)-r  ( I i ,N-1  ) ) 

OF ( I I , l ) = H I* t M II  ,2  )-F { l I , 1 ) ) 

M = N-  l 

rr  i 1 = 2, Nil 

1 r.  F < I1,I)  = .5*HI’<F(II,IM)-MII,!-1I) 
D t TURN 
F NO 


i 


' ft  AN  IV  b lcVEl  21 


SR 


1 (?4< 

DATE  = 74266 


17/33/18 


' 1 
12 

13 

14 

15 

16 
17 
;fl 
:9 

0 

1 

2 

3 

4 

5 

6 


S UR  ROUT  INE  SR ( AP , YH , DM S , SF , RW > 

DIMENSION  SF ( 5 ) * RM  < 5 1 
H I = . 2 

DO  1 1=1, NMS 
PP2=-l . ♦ H 1 

ANF=0. +4.«FSK(0P2,YH. 1,11+0. 
A\M=0.+4.*FSR(RP2,YH, I ,21*0. 

DO  2 J = 3 , 9 , 2 
BP1 =- 1 . +J*H  I 
BP2=-1.+(J+1) *H I 

4MF=ANF+2.*FSR(BPl.YH,I,l)+4.fcFSR(BP2,YH,I,l) 
2 ANP=ANM  + 2. • F S H ( BP1 ,YH,  I ,2)  +4.-FSR ( BP2  » YH, 1,2) 
SF  C I ) =AR+HI*ANF+.  5236 
l RM(  I l--  AR«Hl*AMM*,261tt 
RETURN 
END 


a < 


17/33/18 


MINE  Ml  N 
IMY.LT. 
IF  ( I I .EO 
I F ( II  . E Q 
I c ( I I . f Q 

r.L  ni  ? 

1 IF( Y.GT. 
I F ( Y . G E . 
GO  TP  2 

3 CONTINUE 
C-  S = S I N ( 1 
I F ( I I . FQ 

2 CGNTINUE 
IF(  IJ.EQ 
I F ( IJ  . EO 
RETURN 
END 


FSP ( Y ,YH,  I I , I J > 

YHIGD  TO  3 
. 2) GO  TO  1 

. 1)GS  = S IN (1.8  708-  (Y-l. I/IYH-l. ! ) 
.3)GS=SIN(4.71<!'«*(Y-1.)  /(YH-U  I ) 

YH  ) GS  = 5 1 N ( 1 . £>'08*Y/YH) 

0. I Gb=- S I N( Y*3 . 1418) 


. 5 708* ( Y+ 1. I / I YH  + 1 . ) I 
. 3)GS  = S IM4.  712*  ( Y+  L.  ) / I YH  + 1.  ) ) 


. 1 »FSR=GS*COS( l 5 708*  Y ) 
.2)FSR=CS*SINI3.1416*Y) 


) i 
12 
) 1 

/A 
)b 
) 6 
7 
.8 
:q 
n 

1 

2 

4 

K 

6 

7 

8 

9 

0 

1 

2 
3 

5 

6 
7 

6 

9 

0 

L 

2 

3 

A 

5 

b 

7 


SUHDPiJT  INF  Ct  NVI N ,H, F ,F ! , A 1 , T l, T2, ANS , l J) 

0 IMfc.’S  ICM  F (3  *40)  »FJ  ( 3 » A 0 1 » AN  S ( :> » AO  I * G ( AO ) 

G(  L)=l. 

E l=FXP(Tl*H( 

F 2 = S X P i T2--F  ) 

IFl=-12/(ri*H» 

IF2=-12/(T2*H) 

00  1 I = 2 * N 

1 G( I (=0. 

on  2 1=2,  I 3 1 

2 r,  m = r,  ( i ) ♦ a i * e i * * < i - 1 ) 

00  3 1=2,  IF2 
11=1-1 

3 G( I >=G( I ) M 1 .-Al  ) *E2  * *( 1 - i > 

ANSI  IJ , 1 )=0. 

on  A 1=2, N 

A ANS(IJ,I)=FI(  I J , I )-F I ( IJ,i)*G(I) 

C CO  FIRST  INTEGRARf.E  STEP 

ANSI  IJ,2)='NS(  I J , 2 ) — ,5’H^I F (IJ,l)«G<2)+FUJf2)*G(m 
C 00  CASE  UF  EVEN  NUMBER  GF  BASE  POINTS 

DO  5 I = A , N , 2 

ANSI IJ,I ) = ANS I I J,l  ( I J,l)*G( l)*F< IJ,2)*G( I - 1 * ) 

NP=  l 

AN1=F(  i ».  ?)*G(  I -l  )«-A.*F<  IJ,  3)*G(  I-2M-FI  IJ,  I )*G(1} 

! F { NP . EG , A . GO  TG  5 

NE=NP-2 

DG  £ J = A , NE 

6 AN  1 = AN 1 *2 . * F ( I J , J ) *GC  I- J *• ! ) *A  .*  F (I  J , J v 1 ) *G  I I - J ) 

5 ANSI IJ , I ) = ANS( IJ,  I l-H*ANl/3. 

r.  no  c asf  of  ooo  number  of  basf  points 

00  8 I = 3 , N u 2 

ANl=FIIJ,ll*Gll'H.*F!!J,2l*GII-lhF(IJ,II»Glll 

I F ( I.FQ.3 IGC  TO  8 

NE=I-2 

DO  9 J = 3 » N F ,2 

9 AN l = AN  1 +2.* F I l J , J 1 *G< f - J ♦ 1 ) +A . *F ( I J , J+ L ) “G ( I - J ) 

8 ANSI  I J,  I )=ANS(  I J , ! )- H*ANl/3  . 

200  FORMAT  {•O'  ,2X,'G(  » , I 2.  , • »=•  , E 1 2 - A ) 

RETURN 

FNO 


t r"  wjima'w  \w  n?  v^r^csn  ~ 


1 ;;7 


< 


SAIL-HULL  INTERACTION  RESPONSE  PROGRAM 
INPUT:  1)  Vehicle  Geometry 

2)  Trailing  Vortex  Sheet  Angle 

3)  Exponential  Approximations  co  Circulation 
Response  of  Sail 

OUTPUT:  1)  Step  by  Step  Forces  and  Moments  on  Hull 

due  *o  Unsteady  Sail  Wake 


PiBiajgiww 


pair 


• /^Cf&gvxitrj-je* ne*ex\*v 


BJWTMisa 


IV  G 1 RFLFASE  2.0 


MAIN 


1i;  8< 

OATF  = 74  05R 


1 5/",  2/ 3 


0 I Mf  NS  I r'N  V°  ( AO,  ? » , X(  40  ) ,R  T(  4 0 ) ,OWR(  40*40, 3) 

01  MENS  ION  OWN  A n,  /,o  ,3)  , UF  ( 3 , 4 0 ) , UM  (3  , AO  ) 

D I MFN  S I ON  OUF  ( , 40 ) » HUM  ( 3,40)  ,CUF(3,40),CIJM(3,40) 

PFACIOO , NVP, NMS , NSD ,NX ,»HS , RSL , AP ,TH, AOW, STFRN 
RFA0200,  ( VP ( I ,t)  ,1=1, NV P ) 

RFA0200,  (VP  Cl,?), F = l , N V P ) 
wf  A C l 50  » C » T 1 ,T?, A 
150  rORMA  T ( 4E 1 2 . 4 ) 

X I 1 ) = 1 . 

P T ( 1 ) = P S L 

CALL  WCT(ROW,STFRN,X,RT,NSO) 

PPIM400,  (X(  I ),PT  { I ) , 1=1,  N SO) 

4  00  FORMAT ( *0'  »2X,  ' X = ' ,E12.4,4X, 'P  = ' , E 1 2 . 4 ) 

00  1 1=1, NVP 
xo  = v? ; i,i) 

Y 0=  V P ( 1,2) 

CAIL  UTRLn(XO,YO,Th,X,AR,RT,NMStNX,NSD,ItnWI,RHS,RSL) 

1 CALL  UBNCF XO.YO ,TH, X ,PT , AP ,NMS »NX ,NSO, I ,DWB, RHS,RSL ) 
00  2 1=1, NVP 

PRIM3  00,  ( OWKI  , J,  1 ) , J=  1 ,NX) 

2 PR INT300,  (OWRII , J, 1 ) , J = l ,NX) 

CALL  FM(QwlfnwR,NX,NVS,VP,NVP,X,RT,UF,UM> 

PR  I NT  5 50 


550  cORMA'r  ( '0*  ,2Xf  'HULL  RFSPONSF  TO  UNIT  STEP  OF  CIRCULATION') 
PRINT5  00,(X(I),UF(1,  I )fIJM{  1,  I ),  1 = 1,  NX) 

H = X ( 2 1 - X ( 1 ) 

CALL  OFPYiNX,Htur,CLF,  1) 

CALL  OFRY(NX,H,IJM,OLMf  1 ) 

pp i nt 6 oo, ( x( i ) ,ouf ( i , i ) ,mjM( i , i ) , i =i ,nx > 

CALI.  C0NVfNX,H,0UF,LF,A,Tl,T2,CUF,  l) 

CALL  CONV  (NX,H,DUM,LM,A,T1  ,T?,CIJM,  l) 

00  5 1 = 1 , NX 

CUF( 1,1) =CUF ( l , I )*3 . 14 1 6*C 
5 CUM( l , I) =CMM( 1 , I )*3 . 1416*C 
PR INT56C 

560  FOPMAT ( *0' ,2X, 'HULL  RFSPONSF  TO  CIRCULATION  OF  SAIL  TRANSIENT* 
PR  TNT5  00, ( X( I ) ,C1 F(  1 , I ) ,CUM( 1 , I ) , [=1 ,NX ) 

100  FORMAT ( 4 I 5 » 6P  5 . 3 ) 

200  FORMAT ( 16F5.3) 

300  FHRMAt'  >0*  ,2X, 1 OF  1 2.4) 

5 0 0 FOPMAT('0,,2X,,X=',F12.4,2X,'CY=',E12.4,2X,,CN=',F;1?.4) 

6 C 0 chpma.-F'O' ,2XV  «X  = * ,Fl2.4,2X,*nCY/DS=',E12.4,?X,,nCN/DS  = ,,F12.4 

STOP 

END 


I t A 


I09< 


1 1 R nj  r 


n A T P = l'<rr>’^ 


I S / A 2 / 3 3 


S'j^Tl  T ! \F  nrt-gn  ( v-  , y , , T , X , P T , AR,NMRtivx,NSC,  I C P . nw  3 , P MS  , P S L ) 

M ^’c  \ R ! f \ X ( ? I ) , ■> ^ ( 7 " ) , OWR  ( AO , ?r  , 3 ) , HI  '.  ) , F ( 6 ) , c ( 6 ) , W ( 3 ) 

/ 1 ( A , M ,r  ) ^ . R'-  ( fl'Mf  -F  ) +f  f P ) 
h = x ! 7> -X  U ) 

7>  — 1 . ) «S  IN  ( T ) 

x : i (x^-i.)^rn^(T) 

' 1 l T = 1 t N m «; 

•"  1 N = l * NX 

I <-  ( N.OT.NROK-U  rr  H 
XN= X ( \ ) - l . 

0 T F = « T { K-  ) 

XH  = X ( N ) 

3H  = FR AO ( XH ) 

A<K=Ak'(tJT^-WH)/(PSL-PHS) 

Vrl=VC-RSL4-APtPTT-ARN 
Yril-APN-OT^t^H 
YHT  = APSIN(YH1/APM*.  6366197 
, CONTINUE 

1 F I b .CT  .NSn  jx  (N)=X  ( N-n  +H 
IP  ( N.r.T  . NSn  ) X N = X N + H 

l P ( y; 1 .Gr . YH1 )C0  TC  2 
UT  = . I*  ( 1 • — YHT  ) 

Y1=YHT4HI 

VI  NT  = Fi|RNO(YHT  ,X)  L ,Y0l  , ZC,  XN,  , YHT  , l FIJBNOI  Y1  , X01  » YOl  . 70, 

l ARN,  YHT,  I ) 4FHRN0U  . ,X0l  ,Y0I  , 7 1 , XN,  ARM,  YHT  , I ) 
on  1 J=2,R,2 
Y1=YFT+J*HI 


Y2= V 1 +H I 

3 VINT  = VINT4-2.<  FUR  Nn(Yl,XGl,YCl, 70,  XN, A RN,YHTv!)t4.*rtlh  NO  (Y2,XGl, 
l , 70 ,XK , APN.YHT , I ) 

OWR ( I C P , N , I )=HI*VINT/3. 

on  rr  i 

2 r>IF  = .ft366l<J7*APSINmi)-YMT 
UP=YHT4-2.^0IF 

IF  (IJP.C.T  .1  . |UP=2  .*<  n IF+YHT  )-  1 . 

01 = (UP-YHT ) *.S 
D2=  ( l . -UP  ) •‘.5 
R < 1 > = .2336192 
R ( 2 ) =- 1 .*R(  1» 

R ( 3 ) = .66  12094 
R (4) =-l.*P (3  ) 
p(5)=. 9324 096 
R(6  )=-l  .*»  (5  » 

W ( 1 |=.067°139 
W( 2 )= . 160761 f 
W ( 7 ) = , 1 7 L324S 
00  A 11=1,6 
A 1 = R ( I I ) 

F I I I ) =7  1 ( A 1 , vht  ,l)P) 

A F ( I I ) = 2 1 ( A 1 , I ) P , 1 . ) 

0 1 = C . 

r.7  = :. 

1 "*  6 11=1,3 


no 


< 


! 1 I.  r^-  ’.c  I 11  \ r '■'ATT  = 7 A D 5 B lO/AP/XB 

X ! - c ( > * ! ! - 1 ) 

< 7 - p ( ? - I I ) 

«1  --f-  (7*1  !-!  ) 

w.>=  F ( ?»  I T } 


'•  1-01  ♦ 

( I ! )(  X 1 

, v )1,YC1,/C, 

X M i ARAJ,  YHT  , ! 

1 ) +F  MB -ID  ( x?  , xn  1 

, y:i 

1 

, XA  , U ■'] 

, VF ' T , ! ) ) 

s 

( I I K ( F'M  -Ml  ( W 1 

, X01.YC 1 ,70, 

XN.ABN,  YHT  , 1 

I ) ^imni  k?  , XO  1 

, Y?1 

1 , A \ t A M , v H T , ! ) ) 
jvs •»  1 1 n*  i i = ( r i *Di  ♦n?*n?  i 

1 C"1'  T i \i  I F 

- T I .'A 

: NT 


•'4  '•  0< 

V 1-  <** 


IV  r,  1 ° £ l FAST  2.0 


I ) T R L n 


74058 


l 5/4?/ 37 


Sim « Cl;T  I NT  UTRLD  (X0  ,Y0,  T ,X  , AR  ,PT  ,NMS  ,NX  ,NSD  , ICP,  DWI  ,RHS  , RSI.) 
DIMENSION  X(20),PTt2C),DWH40,20,3),R(6),E(6),F(6),W(3) 

71  (A,B,C)=.5*(A*(f-E)+C+R) 

H=  X ( 2 I - X ( 1) 

7 3=  ( X C - 1 . ) * S I N ( T ) 

X 0 1 = ( X'3-l  . ) *CDS  ( T ) 

on  i i = i , n m s 

DO  1 N = 2 * N X 
I F ( N.GT.NSDJGP  TO  6 
XN=  X ( N ) - I . 

XC=X(N)-.5*H 
XT  = XN-H 

RTC=  (RT  IM  *RT  CN-1 ) ) * .5 
RHC  =HP  AO ( XC  I 

ARN  = AR* ( (RTC-RHC)/(RSL-PHS! ) 

Y01=Y3-KSl +AR+RTC-ARN 
VH 1 =APN-P  TC  + RHC 
YHT=ARSIN( YHl/ARN* *.6366197 
6 CONTINUE 

IF(N.GT.NS0)X<N)=X<N-1)«-H 
IF (N.GT .NSO) XN=XN*H 
1F(N.GT.NS0)XT=XT+H 
iriYOl.GF.YHl ) GO  TO  2 
H I = « l * ( l .-YHT  ) 

Y 1 = YHT+H l 

VINT  = FUTPLD( YHT  ,X01, YOl ,7  3,XN,XT, ARN, YHT, I ) +4  . *FUT  P.L  0 ( Y 1 ,X01  ,Y01  , Z 
1 0,XN,XT ,6RN»YHT,  IMFUTRLPI l , ,X01 ,Y01 , ZO,XN, XT , ARN, YHT, I ) 

DO  3 J = 2 ,8,2 

Yt=YHT+J*HI 

Y2-YUHI 

3 VI  NT  = VI  NT«-2.*fijjrLo  ( Yl  ,X01  ,Y01 , 7-),XN,XT  , ARN,  YHT  , I ) +4  . * FI  IT  P L D ( Y 2 , X 0 
1 l , Y01  , ZO, XN, XT , ARN, YHT,  I I 

PHI  ( ICP.N,  I )=HI*VINT/3. 

GO  TO  l 

2 DIF  = . 63661  97 *ARSIN(Y01) -YHT 
UR-  YHT*-?  .*01  F 

IF{ UP .GT. I . I UP  = 2.* I DIF* YHT )- I . 

Dl=(UP-YhT)*.5 
02=(  I --IIP)*. 5 
R ( 1 ) = . 2386  19? 

P(2)=-l.*R(l) 

R( 3) =.6612394 

R 14 l=-l .*P  <3  ) 

P (5)=. 9324695 
R ( 6 ) =- I . *R ( 5 ) 

H (1) = .4679139 
W(  2 )=.  36076  16 
W (3  I = . 17  13245 
DO  4 11=1,6 

A 1 = P ( I I ) 

E II  I)  = Z1  ( Al , YHT, UP  J 

4 F { I I)  =ZI  ( Al  , I IP  , l . ) 

G 1 = C . 


IV  r,  1 p M C a S F •-> . c 


u t p in 


DAFT  = 7*053 


G?=c. 

fo  5 11  = 1,3 

X1  = E(  2M  I - 1 ) 

X2=  E ( ?* I I 1 
wi=F(?fii-n 
w?=f  t ?*  m 

G L = G 1 + W ( II  ) *(R.ITRLn  ( X l , XO  1 , YO  1 ,Z0,  XN:  XT  , ARN  , YHT,  I)  +RlTRLO(  X2  ,XOl  ,Y 
131 , ZC ,XN,XT, ARM, YH1 , I ) ) 

5 G2=G?+W(I!)*(rUTRLr:m,X0l,Y31.7C,XN,XT,ARNfYHT,I)mjTRLC(W?tX0l,Y 

101 , ZO , XN, XT , ARM , YHT , I) ) 

r>wi  ( i cp,n,  i ) = (Gi*nnr,2*D2 ) 

1 CONTINUF 
RFTURN 
F NO 


11^ 


IV  G 1 Pf  LEASE  2.0 


FUTRLP 


DATF  = 7 A 05  8 


1 5/A2/33 


FUNCTION  FIJTRL0(Y,XC,V0,70,XN,XT,  AR,YH,  II  ) 
IF{n.F0.2)GC  TO  I 

I F ( II .EQ. 1 )GS=COSI l .5708*1 Y- 1 . )/( YH- 1. » )/< YH- L. ) 

IF ( I I .F0.3)GS=3.*CCS (A.  7 12A*< Y- I . )/{ YH-L . ))/  I YH- l . ) 
GO  TO  2 

L IF(Y.CT.YHIGS=CnS(1.5708*Y/YH)/YH 
IF ( Y.GE .0. )GS=-2 .+CCS (Y*3. I* 161 
2 CONTINUE 

ETA=AP*S IN(Y*l .5708 J-YO 
OEN=FTA*FTA+ZO*ZO 

FI  = (XN-XO)/SOKT(  ! XN- X 0 I * ( XN-XOI+DEN) 

F 2 = (XT-XO) /SORT ( {XT-XO ) * ( XT -XO ) +PEN ) 

niTR|_D=-.l?5*GS*ETA*(FI-F2)/DEN 

RETURN 

END 


‘ 


r-’--*.  -I  - . -r  vr  •■ ...  Tr~»T miSSSSZfft 


1I5< 

IV  01  Pp|  FAFF  ?.r  FM  OAT  F = 743S3  lE/<.?/33 

S' IP  R Clll  I ME  FM  ( DW  I , GWB  , NSP>,NWS  , VP,NVP  , X,  RT  ,IJF  , IJM  ) 

D I M E N S I 0 N OW  I (40,20, R » , GWR (40, 20, 3 ) , VP<  40,2 ) , X( 20) ,RT ( 2D  I , OF ( 3, ?C) 

1,IIM{  ?,?C) 

ni  MFNS  MM  AF  (70  ) , AM  ( 20  ) 

HI =VP(? , I ) -VP( 1,1) 

"1°  1 I “ 1 * NVP 
X A = V P ( I ,1  ) 

AF(  I ) = A S L P ( X A ) 

1 AMI  I | = AF ( I ) *XA 
nn  3 k=l,nms 

f F IND  T = 0 + RESPCNSF 

S 1 )w  F = 0w  R ( 1,1  , K ) * AF  ( 1)4--  .*  DWP  < 2,  1,K)*AF(?)4-DWH(NVP,1,K|*AF(NVP) 
St)MM=DWR<  1 , l , K ) *AM<  1 ) 4-4  .<TjWH  i ?,  1 , K )*AM(  ? ) + HWB  (NVP  , l,K)*AM(NVP) 
NFND=NVP-2 
GO  2 I = 3 , N ENP , 2 

SUMF=SUMF  + 2 .*DWR(I  , 1 ,K)  *AF  (I)+4.*0WB  ( I4H  , 1 , K ) *AF((  4-1  1 

2 S'IMP  = SUMM4-?.*nwR(l  . 1 ,K) *AM( I )4-4.*nwR(H-l,l,K)*AP(I  + l) 

U P ( K , 1 )=SUMF*HJ/3. 

UM(  K , 1 ) =SIJMM*HI  / 3 . 

C FIND  RESPONSE  W/TRAUFRS 

On  3 1=2,  NSD 

S'JMF  = SUMFMOWB(  1 , I ,*  l-OWRI  1 , I -l  ,K  1+DWI  ( 1 , l,K)|'»4F(  1 ) 
SUMF=SUMF*4.*(DWB(  2*  t»S  J-DWB  l ?,  1-  1 ,K)  fOWl  ( 2,I,K))*AF(2» 

SIJMF  = SUMF+  (DWfUNVP,  I , K | -PWPINVP,  I - l » K I ♦OW  I (NVP,  l,K))*AF(NVP| 

SUM  M=  S UMM  4-  ( 0W8  ( 1 , I , K l-OWR  U,l-l,Kl*nwi(l,!,K>)*AMm 
SUMM= SUMM+4 . * ( DWP ( 2»  I ,x  1-DWP  (2,J-l,K)*OWI<2,!,K))*AM(2> 

SUE M=SUMMF (DWB (NVP, I ,K I -PW  (NVP, 1-1 ,K » *DW | (NV°, I,K))*AM(NVP) 

DO  4 J = 3 ,NEND,  2 

SUMF=SUMF+2  .*  (DWP(  J , I ,K  1-OViB*,  J,  I l,K'l*DVM  ( J , ,K  ) ) *AF  ( J) 

SUMF  = StJMF*4.M0WP(  J4-1  , [ ,K)-GWB{  JM  , I , K ) 4-DW  I { J «•  l , I , K )>  *AF  ( J + 1 ) 
SUMN=SUMM  + 2.*(DWF( J , I , K I -DwB ( J , I - 1 , K ) 4-QW I ( J , I , K ) ) * AM ( j ) 

4 SUMM  = SIJMM+4.  *(DWR  ( J + 1 , I ,K  |-OWB(  J+l  , I , K)  4-DW  I (Jfl,I,K))*AM(J+l) 

UF  (K,  I ) = SUMF*HI/3. 

3 UM(K,  I )=SUMM*HI/3. 

RETURN 

END 


IV  G 1 KFLfcASf  2.0 


HR  AD 


DATE  = 74058 


15/42/33 


FUNCTITN  HR AD ( X ) 

B0W=-5. 125 
STEPN=7.0625 
X l =-4 . 25 
X2=  .875 
X3  = 6 . 

Rl=.4351 
R2=  .8438 
R 3 = . 3 ? 9 1 

IF ( X.LE.X1 )GO  TO  1 

IFIX.GE.X31G0  TO  2 

A=( P1+R3-2.*P2 )*.5/ (X2-X1 )**2 

B=(R1-R3)/(XI-X31 

C = R 2 

XD  = X-  X2 

P = A’»XD*XD«-B*XD4-C 
GC  TC  3 

1 R=R  1*SQRT( ( X-BOW)/( Xl-BOW)  ) 

GO  TO  3 

2 R=P3*SQRT ( (X-STEPN)/ (X3-5TFRN)  ) 

3 HR AD=P 
RETURN 
END 


A * 

1 i.  f 


IV  G 1 I-  ' L F ft  ^ F J . C 


OATF  = 74  05  8 


1 5/4?/ 38 


s'IHPCIIT  I NF  Wf  T ( A , R , > , P , M 
rl  ! VF  NS  I nN  X ( 70)  , o ( ?C) 

F I NT  { F , F i G » I * H ) = F / ( ( F-  1 .-  ( I - 1 ) *D  ) * ( F-  1 .-(  I - M *n)  «-G*G  ) **  1 . 5 
H=.C?c;*(  P-A) 

HT= ( 0-  1 . ) /N 
N 1 = A ♦ 1 
DO  1 J=1 , N1 
R 0=  « ( J ) 

F 1 = A S L P ( A ) 

A l = A+H 
F 2 = ASL  P ( A 1 ) 

F3=ASLP(P) 

X I N T = F I N T ( E l *A,RO*J*HT)  4-4  . * F INT  (F?  ,A1  »RQ,J,HT  ) + F INI  ( E3,R,R0»J»HT) 
CO  2 I J = 2,  3 8, 2 
A 1 =A+IJ*h 
A 2 = A 1 +H 
F 1= ASL  P( A 1 ) 

F2=ASIP(A2> 

2 XINT=/INT*2.*FlNT(FlfALtPO,J,Hn+4.*FlNT(F2,A2,R0,J,HT) 
XfNT=h*HT*X[NT/L8.85 
R( Jf 1 1 =R ( Jl* ( 1 .♦XIAT ) 

1 X(J*1)=X( JI+HT 
RFTLRN 
ENO 


. ->  rifrS’tT^iViV'1  i ifCtf 


* *^a3B*8sgg^^ 


lib- 


[V  G 1 RFLFASE  ^.0 


4SLP 


DA T F = 74058 


FUNCTION  ASLPIXI 
R0H=-5. 125 
STF  RN  = 7. C625 
X 1 = - 4 . 2 5 
X2  = . 8 75 
XB=f  . 

R L = .4351 
P2=  . 8438 
R3=  .3281 

I  FIX. I F.Xl  )G0  TO  1 

I F ( X.GF.X3)GG  TO  2 

A=( Pl*R3-2.*R?)*.5/<  X2-X1)**2 

R = ( R 1 - R 3 ! / IX  l — X 3 J 

C - R 2 

X0=  X- X 2 

ASL  P = 6. 2832*  (2.  * 4*  A*  XD*  XD*  XD  + 3 . * A *R  *XD*  XTH-  { B*B  + 2 . * A*C  ) *XD«-C*  \ ) 
GO  TO  3 

1 ASLP=3.14l6*Rl*Rl/( X1-B0M) 

GO  TO  3 

2 ASLP=3. L416*R3*R3/( X3-STFRN) 

3 CONTINUE 
RF  TURN 
END 


3 19< 


IV  r 1 p I | r A s F ? . 0 


= 74QSR 


1 S/4?M7 


SU^RHUT I Nf  OFRYIN.H.F.np,  ! I ) 

TI ^FAS I ■>  r M ,40) , CF  { 1,  40  > 

HI  = 1 . /H 

OF(  I I ,N  ) =HT* (F  ( I I f N )-F ( II ,N-  1 ) ) 

OF  ( I I ,1  )=HI*  (F  ( I I ,?  )-F ( I I , 1)  ) 

N 1 = N — 1 

no  1 I - ? * n i 

1 OF(II,n  = .5*HI*(FUI,Ul)-F(II,I-l|) 
RETURN 
E NO 


/■V.--  a, 


" *<Wyv\ 


120- 


IV  C 1 P F l fc  A S F 2 . 0 


DATE  = 74058 


1^/42/33 


SUBROUTINE  CGNV  (N,H,F  , FI,Al,TL,12,ANS,IJ) 

DIMENSION  F (3  ,40) , F I ( 3,40) , ANSI  3,40? ,GI40 ) 

G ( 1 ) = l . 

F l = E X P I T l*H) 

F2  = E X P ( T2*H) 

IF  L = - L 2/ I T l*H) 

I F 2 = - l 2 / <T2*H) 

DO  1 l = 2 » N 

r,  ( i » = o . 

DO  2 I =2 , I FI 
G ( I )=G(  I »«-Al*El**(  1-1  ) 

DP  3 1=2, I E2 
11=1-1 

G ( I ) = G ( I )♦(  l.-Al)*E2+*(  I - 1 > 

ANS I I J, 11=0. 

DO  4 I =2 ,N 

ans  ( i j.  n-F  m j,  i > - f 1 1 1 j,  n*Gi  i > 

DO  first  integraelf  step 

ANSI  I J,2)=ANSI I J,2)-.5*H*(Fl  I J,1 )*G( 2) *F( IJ,2 )*G( L) ) 

CO  CASE  OF  FVEN  NUMPFR  OF  BASE  POINTS 

00  5 I =4  » N ♦ 2 

ANSI  I J , I ) = ANS(  I J * I ) - . S*H* ( F ( I J , 1) *G (I) *F ( I J , 2 ) *0 II - 1 H 
NP=  I 

AN1  =F  ( I J ,?  )*G<  1-1  )*4.*F  (J  J,3)*G<  1-2)  *F  I IJ  , I I^GU) 

1 F ( NP . F 0 . 4 ) GO  TO  5 
N E = N P - 2 

on  6 J = 4 , N F , 2 

AN1=AN1*2.*F ( I J, J»*C< I-J* l)+4.*F ( I J,  J+1)*G(  I - J > 

ANS  ( I J,  n = ANSI  IJ  , I )-H*ANl/3. 

DO  CASE  OF  ODD  NLMBER  OF  BASF  POINTS 
DO  8 1=3, N, 2 

ANl=F(  I J, 1I*G< I )**.*F( I J,2)*Gt I-l  H-F ( IJ, II*Gl 1 1 
I'M  I .F 0.3) GO  T0  F 
NE-- 1-2 

DO  <3  J = 3 , NE  , 2 

ANl  = ANl«-2.*F(IJ»J)*Gti-J-*l)*4.*FIIJ,J  + l)*G(I-JI 
ANSI  I J, l ) = ANS( IJ, I )-H*ANl/3. 

FORMAT  I *0' , 2X,'G( ' , 12,'  J = ' ,F1  2.*  ) 

RETURN 

END 


